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^H ■ Abstract 

"<) . Inspired by the work of Feynman, Deutsch, We formally propose the 

theory of physical computability and accordingly, the physical complexity 
' (-H ' ^ theory. To achieve this, a framework that can evaluate almost all forms 

O ,1 of computation using various physical mechanisms is discussed. Here, we 

focus on using it to review the theory of Quantum Computation. As 
a preliminary study on more general problems, some examples of other 
physical mechanism are also given in this paper. 
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1 Introduction 

1.1 Quantum Computation 

The research of quantum computation has been lasting for about 30 years since 
R.Feynman proposed the concept of so-called 'quantum computer' in 1982 [2]. 
Founding out that there do exist some quantum systems which are suspected 
cannot be efficiently simulated by classical computers, early researchers natu- 
rally speculated that quantum mechanism itself may provide stunning power of 
computation. In order to strictly describe 'quantum computation', researchers 
introduced various new computational models, including Quantum Turing Ma- 
chine J^ and Quantum Circuit Model [5 . However, at that time, no convincing 
evidence was discoverd to support the conjecture that quantum mechanism can 
really be used to speed up the computation of some hard problems greatly. 

D.Deutsch found the first evidence that quantum computers may surpass 
the Turing Machine 9^ in query. In fact, he constructed a special scene in which 
DTM has to query the oracle for 0(2") times to find the correct answer in worst 
cases while QTM need just once query in all cases. 

One of the most remarkable results is quantum factorization, which is due to 
Peter Shor[TUl|TT]. The best classical algorithm for factorization so far has to run 
for 0{log'^ L) steps. However, Shor showed that one can use a family of quan- 
tum circuit, which contain 0(log L) gates and needs ©(L'^ log i log log i) op- 
erations, to get the right answer. 



Grover's Algorithm 'JT. is another successful example of quantum algo- 
rithms. This algorithm can be used to search a database without structure. 
It is easy to proof that the time complexity of this problem for Turing Machine 
is 0{n). However, there does exist a quantum algorithm whose time complex- 
ity is 0{y^). Since it has been proved that this is the optimal algorithm for 
all algorithms that considering quantum mechanics |18| , so the complexity of 
Grover's algorithm can be looked as the quantum computational complexity of 
this problem. 

One of the most important reason that why Quantum algorithms(especially 
Shor's algorithm) seem interesting to computer scientists is that they form a 
challenge to strong Church- Turing thesis, which states that any model of com- 
putation can be simulated on a PTM with at most a polynomial slow down. 

1.2 Physical Computation 

On the other hand, with the exciting research in quantum computation as 
well as other new paradigms of computation(e.g.DNA computation), the idea 
that we may just look physical processes as computations(not just the Tur- 
ing Machine)was also developed. The seeds of this idea can be traced back to 
Feynman[2 , Deuthsch[3j and PitowskyjTj et al. 

It is not very hard to understand and appreciate this idea, for at first glance, 
this point of view has at least three benefits: 

• It can include the concept of algorithms depicted by Turing Machines eas- 
ily, for an algorithm on Turing Machine(its physical implementation) can 
also be thought as a family of physical processes and the corresponding 
measurement. 

• We can try to solve some special kinds of problem with less time or space 
than the lower bound(proofed by one who only consider Turing Machines) 
of them. 

• Being the ones which could be directly simulated, some physical methods 
can also enlighten us to design smart algorithms on Turing Machine. 

What's more, currently, it seems that we cannot exclude the possibility that 
there does exist a family of physical process which can help us to calculate some 
problem which cannot be solved on a Turing Machine in principle. 

However, because of vagueness and extraordinary generality, the theory of 
so-called 'physical computation' has a significent defeat. 

• In many cases, people cannot decide how to define the resource for a 
physical algorithm and as a result they cannot proof or even conjecture 
formally whether a physical algorithm is really superior to any algorithms 
on TM with the same aims. 

Note that the theory of quantum computation is almost free from such defeat, 
for researchers have completed the formal definition of the computational model 
of quantum computation in the early years. Roughly speaking, things tend to 
go wrong when: 

• people adopt a design which the physical axioms it depends on is just an 
empirical one. 



• more than one different systems of plrysical axioms are used. 

1.3 The structure of the article 

Tfie structure of tliis paper is as follows. In Sec. II we introduce two models 
of quantum computation and the definition of complexity respectively. We 
will formally establish the theoretical foundation of physical computation and 
propose the theory of physical computability in Sec. III. In the beginning of Sec. 
IVwe try to use the theory of physical computation to reanalysis the quantum 
algorithms. In the end of Sec. IV, we focus on the topic about how to construct 
problems which take advantage of quantum simulations. 

2 Models of Quantum Computation 

2.1 Quantum Turing Machine 

Quatnum Turing Machine was first introduced by BeniofFT in 1980 and was 
developed by Deustch and Yao. The modern definitioin was given by Bernstein 
and Vazirani in 1997[1]. 

Definition 2.1 (Quantum Turing Machine, Bernstein 1997) Let C* be a set 
of complex nmber a satisfying: For each a, there exists a polynomial time 
algorithm to compute the value of Im{a) and Re{a) close to 2~" within the 
true value. 

A Quantum Turing Machine M is defined as the triple (E, Q, 5), where S is 
a fintie alphabet with an identified symbol ^ , Q is a, finite set of states with 
an identified initial state qo and final state qj ^ qo', 6, the quantum transform 
function 6 : Q x Y. ~> C'SxQx{L,fl}^ rj.^^ qrp^ ^^^^ ^ two-way infinite tape of 
cells indexed by Z, and a single read/write tape head that moves along the 
tape. We define configurations initial configurations and final configurations 
exactly as for DTMs. Let 5" be the inner-product space of finite complex linear 
combinations of configurations of M with the Euclidian norm. We call each 
element phi £ S a. superposition of M. The QTM M defines a linear operator 
Um ■ S ^ S called the time evolution operator of M as follows: If M starts in 
configurations c with current state p and scanned symbol a. The after one step 
M will be in superposition of configurations ip — ^^ atCi, where each non-zero 
tti corresponds to a S(j),a,T^q,d), and Ci is the new configuration that results 
from applying this transiton to c. Extending thi map to the entire space S 
through linearity gives the linear time evolution operator Um- 

Definition 2.2 If Um can keep Euclidian norm, then we say M is well de- 
formed. 

Theorem 2.3 If QTM is in the superposition tp = ^ • aiCi and is observed, the 
probability of the observer gets the configuration Ci is jofip, and then AI is in 
the state -0' = Q- 

Theorem 2.4 We say a QTM is well-deformed if its time evolution operator is 
unitary. 



In QTM, the number of the read/write tape head moves during a computation 
is the cost of time. 

Theorem 2.5 There exists a universal QTM, which is polynomially equivalent 
to any QTMs. 

2.2 Quantum Circuit Model 

The first quantum circuit model was due to Deutsch. Then quantum circuit 
model was improved by Yao[S], who also proved that for any QTM, there exists 
a uniform family of quantum circuit which is polynomially equivalent to that 
QTM. 

Not like QTM, quantum circuit model tends to descript an algorithm by us- 
ing universal quantum gates and circuits without loops. Quantum circuit model 
does not need infinite many quantum gates, but finite many quantum gates 
which called the universal quantum gates. It has been proved that Hadamard 
Gate, phase gate, C-NOT Gate and tt/S Gate are universl. For any finite di- 
mensional U operators, we can always approach it effectively by means of a 
universal family of circuits ^ , which only consists 4 gates above, i.e. 

Ve(3n G '^),EiU,Un) = ma.x\\{U - UnM\\ < e 

W 

The scale of a quantum circuit is defined as the number of the universal 
gates and the depth is defined as the longest path from input to output, if the 
gates is looked as a vertex. 

Both Quantum circuit model and QTM are important models of quantum 
computation. But we do not know whether they are the most natural models of 
quantum computation or do they fully take the advance of quantum mechanics, 
no matter in the theory of quantum computibility and quantum computational 
complextiy. 

3 The Theory of Physical ComputabiUty 

3.1 Observer 

Measurement is in terms of the observer. Though there are many differences 
exist among people's opinions about the exact definition of human beings, we 
prudently assume that an observer is classical, that is, the observer will never 
get incompatible results during one measurement. 

In this article, we will never use the terminology such like 'the observer of 
observer', or in other words, by 'observer' we always mean the last one outside 
the whole experiment. 

In order to unify various forms of results, we require that the observer only 
accept the symbols on a tape (just something like the one of TM) and also only 
use this to initialize an experiment. 

So we define the legal inputs and outputs as the elements in set S^, where 

S = {0,I,*,.} 

and E+ the finite string composed by elements in S. 
The concept of observer is fundamental. 



3.2 Physical States 

We use (usually finite) distinguishable attributes which contribute to the com- 
putation to label the physical states. In addition, though may not be actually 
concerned in computation, three fundamental quantities, namely, space, energy 
and mass are always included in a state for the sake of analysis about resource 
and complexity. 
So we have: 

n C {xil^i X {X2}'^' X • • • X {xn}-^- X {m}®= X {s}® X {e}*^ 

Or more generally(Quantum) 

n C {xil^i X {xa}-^^ X • • • X {xn}"^- X {C"}®' X {C}® x {C}'^ 

For simplicity, fundamental attributes are usually omitted, i.e. 

n C {xj-^i X {X2}^' X • • • X {x„}^" 

For a certain attribute Ai, what really matters is its type which is decided by 
its dimension. Note that dimensionless quantity(e.g. friction coefficient) can also 
be assigned to a null type. When a quantity is expressed by other quantities' 
combination, it's dimension type should be preserved in all cases, or rather, any 
equation should be dimensional balanced. 

For example: 

^[DML'^T-''] .._ ^lD:M]g[D:LT-^]f^lD:L] ^ mgh^^''^''^^'^'^^ 
^[D:ML^T-^] .._ l(-^-)[_D:M](-^2>)[Z):L^T-2] ^ }_(^jj^y2-^lD:ML^T-^] 

are all dimensional balanced. 

3.3 Physical process and their operator o 

Physical process on a state space fi is a state set whose elements are labeled by 
moment t{t G [0, T], T e M+). 

Pee{ (T,P)| reM+,P: [0,T] ^ 17 } 
If two physical processes on Vt satisfies 

(^2Pl)(7riPl) = (^2P2)(0) 

we can define operator o : P x P ^ P i.e. 

^2 o Pi - Pa 
satisfies: 

1. 7riP3=7riPi+7riP2 

2. if < t < TTlPl (7r2P3)(t) = (^2Pi)(t) 

3. if ^iPl < t < TTlPl + 7riP2 (^2^3) W = {TT2P2){t - TTlPl) 

We introduce 1>P<I as the initial state of P, and <lPl> the final state of P, i.e. 
C>P< = {tt2P){0), <P\> = (7r2P)(7riP) 



3.4 Physical Operator and the operation of operator 

Physical operator is the tuple whose first component is a state x in fl and the 
second component is a physical process whose initial state is x, i.e. 

= {{x,P)\xe n,>P<i^x} 

In particular, a deterministic physical operator O means:0 is a physical opera- 
tor, and 

ii t:20{xi) ^ 7r2 0(2:2), then xi =^ X2 

if we only care about the effect the operator do to the initial state, we can look 
operator as a mapping in 17,1. e.O : fi — > O. 

The operation between two deterministic physical operator is defined as 
follow: 

02oOi=02{<lOi{x)t>)oOi{x) 

In more general cases, it is necessary for us to talk about non-deterministic 
physical operators or random physical operators. A random physical operator 
O contains the tuples which has the same initial states but different physical 
processes, i.e. 

6 = {{x,P)\x e n,t>P<:i ^ x}. 

People cannot decide the output 0{x) just by the initial state x € fl. 

Similarly, if just care about cxtensionality, we can look operator as a rela- 
tionship on il i.e.O : fl x il 

We can also define operations between two non-deterministic operators, if 
some preconditions are satisfied. To do so, we first expand the definition of 
some symbols. 

0(2;) = {P\t>P<i^x} 

0{X) = {P\c>P<ieX cn} 

<iO{x)> = {y\y e 0,3P G 0{x)s.t. <] P[> = j/} 

So O2 o Oi (if they are productive) can be defined as 

O2 o Oi{x) = {P2 o Pi|Pi e Oi(x), P2 e 02{<iOi{x)>)} 

Note that all of the random physical operators are exposed to outside world 
by default. 

3.5 Physical Computability 

3.5.1 Deterministic Physical Computation 

Definition 3.1 (Deterministic Physical System) Deterministic Physical System £P is 
a Five- Tuple 



where: 



S = {0, 1, *, .}S+ is the collection of finite string formed by elements in S , 
r2s+ is the the set of physical implementation of E+. 



• fl ~ {ipi,i G A} J7 7^ rig+A is an index set fi is a set of distinguishable 
physical statcs(labeled by their attributes). 

• V : ilj2+ -^ reinitialization operator 

• H : il ^ ft evolution operator 

9 A : il —> 172+ Measurement operator 

Since Hilbert's 6th problem has not been resolved yet, i.e. the whole theory 
of physics has not been axiomatized, we do not know that whether there exist 
some additional fundamental mathematical constraints should be included in 
this theory. Now, maybe the only restricitions here are the finiteness of the 
resource cost by a physicall process and the finiteness of the attributes used to 
label a physical state set. 

As a result this system may looks looser than many classical computational 
models and may contains the ability to surpass all these models. We would like 
to let physicists to add more necessary restrictions into the system. 

Of course, when it comes to a specific branch of the physics, we can always 
know what is a legal states and processes. However, we wish to keep some 
freedom, i.e. to let the observer combine various axioms in physics so as to 
optimize the computations. 

Definition 3.2 (Partial Physical Computable Arithmetic Functions)For any 
partial arithmetic function / : N — >■ N is said to be partial physical computableii 
and only if there exists a Deterministic Physical System 

^ = (f7,E,V,H,A) 

which satisfies 

If a; G dom{f), then, 

(Ao-HoV)(a;) = /(x) 

Similarly, we can define Total Physically Computable Arithmetic Functions 

Definition 3.3 (Total Physically Computable Arithmetic Functions) For any 
total arithmetic function / : N — > N is said to be Total Physically Computable, 
if and only if there exists a Deterministic Physical System (il, E, V, H, A) which 
satisfies 

Vx G Nwe have: 

(AoHoV)(.t) = /(x) 

In order to extend the definition of physical computability to non-arithmetic 
functions, we should take into consideration the precision of the measurement 
and computation. Therefore, we need a distance function to measure the pre- 
cision of two values and define the computability as the ability of computing in 
any desired precision. 

Definition 3.4 (Partial Physically Computable Functions) Given a partial func- 
tion / : A ^ Band a metric T> : B x B ^^ M/ is said to be partial physically 
computable with respect to the metric T) , if and only if for any e > there 
exists (17, S, V, "H, A)s.t. for any x G Awe have 



if X G dom{f), 

2?((Ao-HoV)(a;),/(a;)) <e 

Similarly, we can also define Total Physical Computable Functions. 

Definition 3.5 (Total Physically Computable Functions) Given a total func- 
tion f : A ~> Band a metric V : B x B ^ M/ is said to be partial physically 
computable with respect to the metric V if and only if for any e > there 
exists (rj, S, V, H, A), such that, 

VxeA 

p((AoHoV)(a;),/(x)) <e 

3.5.2 Non-deterministic Physical Computatation 

On the other hand, many physical processes are considered to be non-deterministic, 
which enable us to implement so-called 'randomized algorithms' and 'quantum 
algorithms'. Our Probabilistic Physical System is defined as follows. 

Definition 3.6 (Probabilistic Physical System) Probabilistic Physical System ^* is 
a five-tuple: 

^* = (f7,E,V,-H*,A) 

. where, 

• S = {0, 1,*,.}. 

• n^{i:„ie A}. 

• V : ilj2+ -^ r^which is also called initialization operator 

• H* : ilxft which also called evolution operator, which is non-deterministic. 

• A : 51 ^ 175]+ J which is also called measurement operator. 

Non-deterministic does not necessarily cause probability, but let's convention 
that in this article we always discussed the randomness which has a probabilistic 
distribution. 

Definition of the computable functions by means of ■^* is an analog to that 
of ^. As an example, we define Total Non-deterministic Physical Computable 
Functions. 

Definition 3.7 (Total Non-deterministic Physical Computable Func- 
tions (Las Vegas)) For any total fucntion / : N ^ N is said to be total 
non-deterministic physical computable function, if and only if there exists a 
five-tuple 

(f7,E,V,H*,A) 

s.t. 

Va;((Ao-H*oV)(x) = /(x)) 



Because of randomness, for any identical inputs x, the system may call 
different process to compute. The above definition is the counterpart of the 
definition of the so called Las Vegas algorithm. 

Definition 3.8 (Total Non-deterministic Physical Computable Func- 
tions (Monte Carlo)) For any total function / : N ^ N is said to be total 
non-deterministic physical computable, if and only if there exists 



s.t. 



where 



(f7,E,V,-H*,A) 

Pr{i<^ueni(x)occurs} > 2/3 
Event{x) = ({A on* o V)(.t) = f{x)] 



3.6 Estimation of the Complexity of Physical Resource 

For the physical systems defined above, we can even ignore that whether there 
exists a physical mechanism in reality to implement it. Any functions which 
could be written as the composition of the three operators would be considered 
as computable(deterministic version). 



n 


-^H 


n 


vt 




a; 


^2+ 


->/ 


^s+ 



But any experiments which implement a certain system will cost resource. 
We will focus on four kinds of resource, namely, time, space, energy and mass. 

Definition 3.9 ( Resource ) The resource of a physical process IH includes: 

• T: The (expectation of the)total time the whole process consumed; 

• &: The maximum of (the expectation of)the space the whole process con- 
sumed; 

• 9H: The maximum of (the expectation of)the mass the whole process con- 
sumed; 

• €: The maximum of (the expectation of)the energy the whole process 
consumed. 

andfR= (T,6,OT,€) 

In the above definitions, the metric of them could be selected as the common 
ones. Today, most physicists tends to believe that mass and energy are not 
independent, neither do time and space. But for convenience, we still focus the 
primitive forms of resource, for actually we don't care about the independence 
here. 

In the above definitions, we don't talk about the potential possibility that 
even time could be reused. 



10 



We convention that the resource is with respect to an inertial system, i.e. 
the observers obtain their results when they are in an inertial system to the 
system running the 'algorithms', so as to rule out the paradoxes because of the 
theory of relativity. 

In many cases, we just cannot get a infinite precise estimation about the 
resource, but for our purpose, we actually do not need such things. Of course, 
there may exist some cases when we could not get an estimation without any 
promise of any precision, however, we will not use such processes to construct 
our implementation. 

The projections of the fundamental attributcs(resource) are t:<xii{S), 7rg(5) 
and 7r£(5) 

The resource a physical process consumed is: 

IP = TTlP 

mP = max{7rOT(S'),5e7r2P(0,7riP)} 
eP = max{7re(S'),5G7r2P(0,7riP)} 
eP = max{7re(S'),5e7r2P(0,7riP)} 

In general cases, when we have to discuss the process of superposition, the 
resource can be defined as: 

TP = TTlP 

mP = max{E[7raR(5')],S'e7r2P(0,7riP)} 
<&P = max{E[7re(S')],S'e7r2P(0,7riP)} 
eP = max{E[7re(S')],S'e7r2P(0,7riP)} 

So it is easy to see that 

TP20P1 = 7riPi+7riP2 

97IP20P1 == max{OTPi,a«P2} 

e P2 o Pi ^ max{ mPx , a7lP2 } 

€P2 o Pi == max{ mPx , 971P2 } 

According to the definition above, the resource of a non-deterministic physical 
operator O which is initialized by x € il should be defined as: 

TO(x) = E[TO,(x)],0,(x) e 0(a;) 

S«0(x) = E[OTO,(x)],Oi(x) € 0(a:) 

<SO(x) = E[eO^{x)],0^{x) eO{x) 

<BO{x) = E[€0,(x)],0,(x) e ©(x) 



So for operators' operation, we have: 



T02oOi(x) = TOi(x)+X02(<lOi(a;)[>) 

m02oOi{x) = max{9JlOi(x),9Jt02(<lOi(x)[>)} 

6C'2oOi(x) = max{60i(x),602(<lOi(a;)[>)} 

e02oOi{x) = max{£C'i(x),€02(<lOi(x)[>)} 
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3.6.1 Framework for the Complexity with respect to General Phys- 
ical Resource 

Definition 3.10 (Resource(deterministic))A resource the physical process 
which complete the whole computation consumed $H^ including: 

1s^ = %i{AoHoV)ix)) 

6^ = 6((Ao-HoV)(a;)) 

Mg^ = 2n((Ao-HoV)(a;)) 

e^ = (£((Ao-HoV)(a;)) 

i.e. d\,^ = (1,32, 6,32, 9H ,3a, €,32) 

Definition 3.11 (Resource(Las Vegas))A resource the physical process which 
complete the whole computation consumed d\^ including: 

T^ = %{{AonoV){x)) 

6^ = 6((AoHoV)(a;)) 

m^ = an((Ao-HoV)(a;)) 

e<32 = €((Ao-HoV)(a;)) 

i.e. *K,^ = (T,^, 6,32,971,32, 2; ,3a) 

Definition 3.12 (Resource(Monte Carlo))A resource the physical process 
which complete the whole computation consumed d\,^ including: 

T,32 = T((Ao-HoV)(a;)) 
&3^ = 6((Ao-HoV)(a;)) 
OT<32 = 2)l((Ao-HoV)(a;)) 

e<5» = (£((Ao-HoV)(a;)) 

i.e. $H,32 = (1,32, 6,32, 9H ,32, 6:32) 

The corresponding concept of complexity should be defined as the resource 
consumed with respect to the length of the input. 

Definition 3.13 Complexity The complexity of a kind of resource is a func- 
tion of the length of the input x n = [log x~\ 

Conplexity^ (n) = inax{T,32(j.)|n — 1 < logo; < n} 

Conplexity^ (n) = max{9Jl3>>(2,)|'^ ^ 1 !i log a; < n} 

Conplexity^ (n) = iaa'x{& ^/x\\n — 1 < logx < n} 

Conplexityg (n) = max{(£^(2,)|n — 1 < logx < n} 
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3.6.2 Some Common Examples 

It is interesting to find some new methods to compute problems without the 
help of universal Turing Machines. Though the ages, people have found a lot of 
such examples, the most famous of them are: 

• Measure the volume of an object by putting it into the water; 

• Obtain the centroid of an object by two suspension method; 

• Compute the function sine by analog circuit; 

• Decide the path of minimum cost using Fermat' Principle; 

• Calculate the mean of numbers by the second Law of Thermodynamics. 
Actually, we can give even more similar examples: 

• By making use of resonance, we can easily find the desired tuning fork 
from a heap of tuning forks. Otherwise, we have to look up the label of 
them one by one and even have to compute the frequency one by one if 
there is no labels on them. 

• We can compute the square roof of an given number x by the law of free 
fall. Prepare a vacuum tube T of length x and let if stand vertically, then 
let an object o which is small enough fall. Get the time t when it touch 
the bottom, and we have ^/x — t/c^ where c — (2/g)^'^. 

• We can sort a series of numbers by dangle poises by strings, where the 
strings satisfies Hooke's law. Given an array of numbers {xi\ construct 
or find poises whose mass is just a;^, then dangle them by strings with the 
same stiffness coefficient. When the system is stable, the position of the 
poises with respect to their weight just indicate the relationship desired. 

However, it is hard for us to estimate the cost of the methods above just 
after we describe them informally. So we select a part of them to analyze next. 
Conventions:^; is the representation of number in digits, [x] is the value of x, 
means the attribute A has the value x. [x\^ means the representation of 



\x 



A 



quantity x though not on the tape. 

3.6.2.1 Mean of Three Numbers 

Given three numbers, compute the mean of them making use of law of thermo- 
dynamics. This idea comes from Pitowsky [7]. 

The strict description of the problem: Given three numbers xi, X2, X3 G 
[ 0, 100 ], compute 

X = (precise to 2 decimal digits). 

Pitowsky suggests that since all of the three numbers less than 100 and bigger 
than zero, note that the freezing point of water is C° and the boiling point 
of water is 100 C° under the one standard air pressure. So for each number Xi, 
we can prepare the corresponding water of volume V and temperature of XiC° . 
And then pour the water of three vessels into a bigger one, whose volume is 
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V'{V' > 3V), and wait. After the water arrived at the balance point, measure 
the temperature. Of course, we assume that during the whole procedure, no 
calory is lose. 

Apparently, the physical state the method above deal with is the temperature 
of water, so we have 

n = {t\tie [0,100],i = l,2,3}^. 

on the other hand, we suppose the water is heat up from C°, i.e. the initial 

state of the experiment is ([0]"^, [0]"^, [0]"^). 

Therefore, the process could be depicted as following: 

V : E+ — > O, heat up the water to the desired temperature 

V(X1,X2,X3) = ([Xi]^, [X2]'^, [Xs]'^) 

"H : O — > O, admixture the water of different temperature, the second law of 
thermodynamics is used 

A : O — > S+, measure the temperature of the water 

Ai[xf,[xf,[xf)^x 

For this problem, since the precision is finite, and there are only constant (three) 
numbers and the numbers are bounded, we can easily deduct that $H,^ is a con- 
stant. As a matter of fact, for Turing Machine, we can also find a constant 
resource costing algorithm which is just looking up a finite list to solve the 
problem. 

3.6.2.2 Sorting Without Repeat 

Description of the Problem: 

Inputs: Finite number series of length n: 

A^ {x^\xi e Z+n [0,M](0<i <n)}; 

Outputs Finite number series of length m: 

B = {xj\xj e ^(0 <j< m)}, 
s.t. if ji < J2 then Xj-^^ < Xj^ . 

Our plan is: for the given series, select a series of poises of length n, s.t. the 
mass of the ith poise is equivalent to the ith number. Dangling the poises from 
right to left by strings, whose restoring coefficient are k. Wait until the system 
is stationary, open the parallel light source and measure the projection onto the 
vertical ruler at the right end. The measurement could be done by machines 
and present the results onto the tape for observer. For Example, we can embed 
some photoconductivc diodes in the ruler by graduations, diodes who is not 
triggered should be read. 
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The physical state the method is primarily concerned with is the mass of 
poise M, the horizontal positions of the poises X and the vertical ones Y, the 
projections Y' and the boole value B indicating which diodes is triggered, i.e. 

So we have 

V : E+ -^ O (Selecting poises) 

V(®r^O = {®^[x^]''[^fm'' ®f [(.7,0)]^' ><^) 
"H : O ^ il (Dangling poises) 

n' -M^ f7(0pen parallel light) 

A : f7 ^ S+(read the projection) 

A{®^{[x,r[jf[[x,]g/kf ef [(j, J e yl)]^'x^) = e,,(x,0 

satisfies if ji < J2 then 

Considering the ideal implementation, we conclude that the Uig^ is linear, 
which is superior to Turing Machines using comparisons, for the complexity for 
them is proofed to be O(nlogn). However, there does exist Turing Machine, 
which is not based on comparisons, also has a linear time cost. 

Note that if the number series is boundless, the complexity of the method 
above will be exponential. This is the common defeat of most analog computers. 

3.6.2.3 Volume of irregular shape 

For this issue, we shall restrict the range of the saying 'irregular' so as to rule out 
the object with infinite length of description. So actually, we tend to discussed 
a subset of the set of all irregular shape. 
Description of the problem: 

Inputs: point series of length n:(xi,yi)(l < i < n), satisfies c + r < Xi < 
a — c — r, c + r < yi < b — c — r 

Outputs: The volume of the box of length a and width b and height /ig, not 
including the scries of cylinders which is induced by the series of points. 

Our plan is simple. Assume we have a box of material of dense p, and a punch 
to extract circles from it. Then we measure the mass of the rest then divide it 
by its dense or just put it into water. V : S+ — > f2 

V(©r=ia(a;.,y.)) = [phoiA - U^o,)]'' [of 

HiiphoiA - U^^,c,)]''[Of' = [pho{A - ii7=iC^)]''[phoiA - Li^^,Ci)f 
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A : O ^ S+ 

A[pho{A - U^^,c,)]^[ho{A - ur=iQ)]^' = ho{A - U^^.c,) 

Apparently the resource complexity for this method is linear with respect to 
the number of the points. However, because most people think that wc cannot 
do infinitely measurement during one experiment, this method can only provide 
the result of finite precision. This is a good news to Turing Machines because 
this implies there exists a Turing Machine which is almost equivalently efficient. 

This may be astonish to someone, who may thought that a TM should at 
least solve the equations first. However, because of the finite precision, Turing 
Machine can just split the object into lattice and use the so-called scan-line 
algorithm to compute the algorithm. 

3.6.2.4 The centroid of Irregular Shape 

Just as the last example, we restrict our topic into the same subsets of all cases. 
Description of Problem: 

Inputs: point series of length n:{xi,yi){l < i < n), satisfies c + r < Xi < 
a — c — r,c + r<yi<b — c — r 

Outputs: The centroid of the box of length a and width b and height ho, not 
including the series of cylinders which is induced by the series of points. 

The method we suggest is similar to the last one, the difference of them is 
that this time we will record some points. 

V : s+ -^ rj 

'^m=Mx,,y,)) ^ [ph,{A~UUc,)r[0f[0f[0f' 

Hib/io(A-Ur=iC.)]^[0]^[0]^'[0]^" = [p/io(A-Ur=iC.)]^^[f^+yo]^[0]^'[0]^" 

|c- Vol 



-H2 : rj ^ rj 



Mr *- ^ I T^lSrnlS'miS" _ r„;. fA 11" „ MMr "^ ^ ii^lSr'' ^0 , T/nS'rnlS" 



H2[pho{A-U-^,c,)r[j--^^+VonO]^ [0]^ = [pho{A-U-^,c,)r[j—^^+Vonj—^^+V^]^ [0] 



. c- 


-Vo 


-\c- 


-K,| 


c - 


-Vo 



[phoiA UlLiQ)]*' [j--^ + Von-—^ + V, 



\c-Vo\ "' '\c-V^\ 



oJ 



c 



A : r? ^ S+ 

4-11" .rA]'^^, , 



AiphoiA ur.icO]^^[^^ + Vof[^^ + V^ficf" = c 
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The time the system cost from oscillating to stillness can be bounded by a 
constant. Because of the same reason this method does not break up the lower 
bound of Turing Machine. But for some other things, we tend to pay more 
attention to it. Some relevant issues will be discussed in Sec-V. 

3.6.3 Graph Isomorphism, Graph Spectrum and Oscillators 

In this part of the section, we shall talk about a complex example in detail. We 
do not mean to show that the method we designed here is superior to all of the 
TMs constructed by the people of the same aim. We just want to demonstrate 
a new style of computation. 

3.6.3.1 Spectrum of Graph 

Suppose X ~ {V, E) is a graph, A is it's adjacent matrix. We say /a (A) is the 
characteristic polynomial of X, also denoted by /x(A). (Ai, . . . , A„), the whole 
root of /(A), is called the spectrum of graph X. 

Actually two different adjacent matrices may represent two isomorphic graphs. 
If we alter the permutation of the number of the vertices, A will become P~^AP, 
where P is the corresponding permutation matrix. However, the characteris- 
tic polynomials of them are the same. Therefore, /x(A) and the spectrum 
spec(X) = (Ai, . . . , A„) are uniquely determined by X. 

For the relationship between spectrum and graph, people conjectured that 
graph can be uniquely determined by spectrum, i.e. suppose 

spec(A) — spcc(i?), 

can we conclude that 

Unfortunately, the different graphs of the same spectrum were found soon. 

Nonetheless, calculating the spectrum is also important. Because we can 
know a lot of crucial properties, such as the extcnsionality, rapid mixing time 
of Markov chains on the graph, by the spectrum of the graph. What's more, 
when two graph have same spectrum, and spectrum is never repeating, we have 
a polynomial time algorithm to check whether they are isomorphic. 

1. Input graphs G1G2, compute their spectrum, denoted by A1A2. 

2. Compare the spectrums, if Ai 7^ A2, then return NOT ISOMORPIIIC;else, 
continue; 

3. Check whether the product of the two similar matrices is a permutation 
matrix, if it is return ISOMORPHIC, otherwise return NOT ISOMOR- 
PHIC; 

Notation: Here by Ai 7^ A2 we mean after sorting their eigenvalue, the two 
series are not identical to each other. And accordingly Gi, G2 should also be 
altered into Gi, G2. But for convenience, we do not differeciate Gi and Gi. 
Proof: 
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If Ai 7^ A2, then Gi '^ G2- So we only eonsider the case in which Ai — 

A2=A. 

i.e. Suppose 

then we have 

thus 



Gi = PAP^, G2 = QAQ^ 
P^GiP = A = Q^G^Q 



Gi^{QP^fG2{QP^) 

by the preconditon,A is never repeating, so PQ is the unique orthganol 
matrices. 

the rest is to show that if Gi, G2 is isomorphic, then QP"^ is the permu- 
tation matrix desired. 
In fact, if Gi = G2,then there exists a permutation matrix S s.t. 

Gi =^ o G20 

Since G2 = QAQ^, the formula above means 

Gi = S^QKQ^S = (g'^5)'^ A (Q^S*) 

Because of the uniqueness of P, we can conclude that Q-^S* = P"^, and by 
orthgonality of Q, we obtain 

s = gp^. 



n 



3.6.3.2 Harmonic Oscilators of multi-freedom 

Suppose s is the number of freedom of the system, qao{oL = 1, 2, . . . , s) is the 
general coordinates when the system is in balance. Without lose of generality, 
we can always assume that q^f) is just zero, i.e. q^f) ~ Q{a = 1, 2, . . . , s). 

Because we only talk about little vibration, so we only keep several terms in 
the Taylor series of the Lagrangians L of the system about qao- 

The potential energy: 

Note that Vq can be omitted. Introduce the notation kafi 



rla'^qfj / Q 

which is called the strength coefficient. According to the formula ( ^^ ) = 0, 
the second order of the potential energy could be represented as 

^ = 2 Xm kapq^qp. 
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Then assume r^ = fiiq) is not relevant to time, i.e. the obhgation is constant, 
so the kinetic energy is: 



T = 



1 " 

IE 



n s s 



rriiri ■ r,, = 



1 ■^-^ ■^-^ ■^-^ wi I lyi I . 



t=l 



i=l a=l (j=l 



la dqp ' 



Introduce the symbol ttIq^, 



Edri dri 
. dqa dqa 

then the kinetic energy could be represented as 



T 



a=l/3=l 



mapqaqfj. 



Keep the formula above to second order and since qaqfi is second order map 
should be expaned to zeroth order. In other words maf) could be looked as 
constants, we just take the value of them when the system is in balanced point. 
So the Lagrangian could be written as 

-^ = 2 E '^{'^afiqaqfi - kapqaqfi)- 

Thus the Lagrangian equation is 



d d /I 
dt dqa \ 2 



s s \ rl I ^ ^ ^ \ 

Yl E "^PiW-r - 7^ ^ 9 E E h-rlpq-r = 0- 
/3=i7=i / '^'^" V %=i7=i / 



di (2E"'"7 

7=1 



q-/ 



1 ^ \ /i ^ 1 

2 Yl ™/3"9/3 + 2 E ^"-r^T + 2 

,3=1 / \ 7 = 1 



13=1 



'f3aqi3 



therefore 



Let 



y^ mafsqi^ + y^ kg/iqp = (a = 1, 2, . . . , s). 
/3=i /3=i 



g^=A^e^* (/3-l,2,...,s). 
Take it into the former formula, we get the hnear equations for Ap. 

s 

y^^jmapX^ + kap)Ai3 = (a = l,2,...,s). 

,3=1 

If the equations have non-trivial solutions, then following conditions should be 

hold: 

TOiiA^ + fcii mi2A^ + fci2 ••• misX'^ + kis 

TO2sA2 + k2s 
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m2lA2 + /C21 '7l22A^ + k 

nisiX^ + fcsi ms2A^ + ks2 



rUssX^ + kg 
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This is the equations of times s of A^, and we can get s A^, denoted by 

Xf G-l,2,...,s). 

3.6.3.3 The characteristic oscillators for a Graph 

Making use of the conclusions above, we construct a specific oscillators for any 
given connected graph. 

Denote the vertices of graph by numbers 1 ^ n, according to any order. 
The mass of a vertex is set Ig. Connect the vertex 1 and n to ends by strings 
whose k is zero by that direction. For the rest, we connect them according to 
the adjacent matrix, i.e. if Aij = l(Note that Aij ~ Aji), connect vertex i 
and j by a string whose k = 1. Let's study the motion of the system: First, 
if two vertices is not connected by string, we have fc^^ = fc^Q, = 0. Second, 
the vibration is little, so string is not an obligation. And we take the general 
coordinates as the usual displacement vectors, so niafj — mpa — Sa/s, where Sa/B 
is the well known Kronecker notation. 

At last, we obtain the determinant as follows, which is the characteristic 
polynomial of our system. 



A^+di -Ai2 ■■■ -Ai 
-A21 X^ + d2 ■■■ -A 



2s 



-Asi -As2 ■■■ A^ + rf, 



= 



It has been proofed that A^ < 0. So let —A = A^ ,we can see that the determi- 
nant above actually compute the spectrum of A' which is converted from A by 
adding multi-loops(the number of degrees). If a vertex is in the characteristic 
position, it will take part in the vibrations of all frequencies, if no one is in the 
characteristic position, then they just vibrate with respective frequency. In both 
cases, we'll measure the frequency and differentiate them by means of FFT, so 
as to get the spectrum of A'. 

Apparently, adding multi-loops is not harmful to the decision of whether A 
and B are isomorphic, for if A' ^ B' , then A ^ B If A i2 _B, then A' ^ B', 
which will also be checked by the oscillating system. 

3.6.4 Steiner Tree Problem 

Steiner Tree Problem is a problem in combinatorics. The general version of 
Steiner Tree Problem is NP-complete, which implies that this problem is unlikely 
be solved in polynomial time. 

This problem is similar to the Minimal Spanning Tree Problem in metric 
space. The difference is that Steiner Tree Problem allow people to add new 
points v'{v' ^ V) and new edges e'(e' ^ E) into the original graph G, if neces- 
sary. When \G\ = 3, the new point (in this case, at most one point is needed)is 
called Fermat point. 

At a time, some people became to believe that the experiments of soup 
membrane can be used to solve the Steiner Tree Problem. In fact, when \G\ is 
small, say, less than 5, this method really works. However, when the number of 
vertices is 10 or more, this experiment just cannot give the right answer. One 
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can attribute the failure to different reasons and derive various explanations, 
among which that it is just the errors in the experiment cause the failure made 
many people conjecture faithfully that classical mechanics can be used to solve 
NP-complete Problems in polynomial time(So they try to proof P=NP). 

In fact, the foundation of the experiment is the well-known property that 
the membrane will be stay in a stationary state, where the surface it produces 
will be just the minimal surface. Unfortunately, this theory has nothing to do 
with the fact that the membrane can arrive at the stationary state fast. What's 
more, no one can proof the soundness of such property under the framework of 
classical mechanics. 

3.6.5 DNA Computation 

In 1994, Adleman used a probabilistic DNA algorithm to solve HP problem (Hamilton 
Path Problem). HP problem is NP-complete, which implies it is difficult to find 
a polynomial algorithm to solve it ll4j . 

In order to understand Adleman's method, the following knowledge seems 
necessary. 

(1) DNA contains chains consisted by four types of nucleotides, denoted by 
A, C, G and T. 

(2) These nucleotides forms complementary couples, i.e. A and T are comple- 
mentary, C and G are complementary. If the corresponding positions of 
two DNA chains are complementary, they will patch up as the twin-helix 
structure. 

(3) PCR, which proposed by Kary MuUis, is method to reproduce the specific 
chain we need. 

(4) There is a machine called 'sequencer' which can be used to read out the 
series of a DNA chain. 

Adleman's Algorithm contains five procedures (Suppose \G\ — n): 

(1) Randomly produce the paths in the Graph, encoded by DNA chains. 

(2) Keep only those paths which begin with Vin and end with Vout- 

(3) Keep only the paths whose length is n 

(4) Keep only those paths which enter all vertices in G at least once. 

(5) If any paths remain, return 'True', else return 'False'. 

Note that the first step of Adleman's Algorithm which is usually thought to 
be work as an initialization operator V is not polynomial with respect to the 
resource mass 9K and space 6 at least. Considering asymptotically we can 
only sequentially get the mass the algorithm need, so actually 0{n\) mass can 
cause 0{n\) time T. As a matter of fact the other steps of this algorithm, 
which require exponentially molecules fully blend by polynomially increasing 
contacting facades, also cost a lot of resource 1. 

It is not very hard to appreciate the conclusion that we can obtaining great 
power of computing suppose we are provided with corresponding quantity of 
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mass, and do not take the cost of preparing such equipment at all. For one 
thing, let's consider the following ideal model. 

Suppose we have enough universal Turing Machines, each of them are de- 
noted by their footnotes. What's more, by some altering in the definition, these 
UTMs have the ability to transmit their results to others. And the condition 
of two UTMs Ui, Uj{i ^ i) could communicate to each other is that they are 
adjacent to each other, denoted by Adj{Ui, Uj). 

So the computational model constructed following, called 'Turing Tree', can 
exponentially speed up the computation of any NP-complete problems. 

Definition 3.14 (Turing Tree) Suppose we have infinite many UTMs, each of 
them denoted by unique footnotes, and 

Adj{U^, Uj)-^j = 2i + l\/i = 2j + l\/j^2i + 2\/i^2j + 2, 

then we call this Turing Tree. 

It is easy to see that the following relation holds: 

Adj{Uo,Ui),Adj{Uo,U2) 
Adj{Ui,U3), Adj{Ui,Ui),Adj{U2, U5),Adj{U2, Ue) 



For example, a TSP problem can be solved as following: 

a The Observer input the weighted complete graph G to the Uq, Uq decode 
to a permutation and compute the sum of the weight, and then transmit 
G and flag i^==0 to Ui,U2. 

h For index iAfter Ui get F — and G, it check whether i < [log2 n!] ,if 
the answer if 'yes' then decode i to a permutation and get the sum, and 
transmit G and F = to U2i+i and [/2i+2;else check whether i — [logj n!] , 
if it is true, decode i to a permutation and get the sum, then submit the 
weight sum to the [/ri/2]-i- Else, do nothing. 

c For index i, after Ui get F — 1 and two sum(come from U2i+i,U2i+2), 
if it's index is not zero, then submit minimi, 6*2^+1, S'2i+2} and i^ = 1 to 
U[i/2']-i- Else return min{5o. Si, S2} and write it on to the tape. 

It is easy to check that the subprocedure of the algorithm which is used to 
decode a natural number to a permutation is polynomial. So the cost of time 
the Turing Tree consumed should be 0{2 log2 n!) < 0(2nlog2 n) (Including once 
sharing the task and once championship for the most value). So it is the time 
to answer how to easily construct a big enough Turing Tree. 

3.7 Preliminary Discussion of the classic theory of Com- 
putation 

3.7.1 Turing computable is physical computable 

The subject about the existence of a theoretical physical system which can 
provide an implementation of universal Turing Machine has been studied by 
many scholars. In addition to the current implementation of computers, scholars 
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have constructed many other wonderful designs on various axiom systems of 
physics(e.g. Classical Mechanics, Quantum Mechanics). 

Of course, the results above only imply that it is the ideal mathematic model 
for a family of physical phenomenons can be look as equivalent to UTM in terms 
of computability. After all, we cannot know for sure that some theory of physics 
is totally right. Because of this, when we talk about the ability of computation 
for a certain family of physical system, we always assume either of the two 
preconditions following: 

• The ideal mathematic model of some branch of physics is believed to be 
absolutely right. 

• At least in a very large scale, the theory works. 

3.7.2 PLATO Machine 

For several decades after the Church- Turing Thesis was proposed, people failed 
to find a counter-example of it. This kind of counter-example, if they really 
exist, should satisfies the property that most people think they can be effectively 
computed in principle, and no Turing machine can compute them. 

However, many physicists tend to make efforts in another direction, that is, 
they want to find a family of processes in nature, whose functional expression 
may not be intuitively computable, nor Turing Computable, but it can actually 
compute a function by measurement. 

Suppose the problem we attempt to deal with now may cost infinite many 
steps for some computational model(e.g. Turing Machine), does it necessarily 
mean that we have to wait infinitely long time to get the results? This is 
not always the case, PLATO Machine, which was proposed by H.Weyl[7], is 
just a counter-example. Though it is named after 'Plato', the designer's main 
inspiration comes from one of Zeno's Paradoxes. 

Specifically, PLATO Machines use (1/2)" seconds to execute the n—th step. 
For instance, suppose the decision problem we want to solve is 3nP(n), where P 
is a predicate and P{x) is used to describe some properties of x. Then PLATO 
machine V will check whether P(l) = 1 holds in 1/2 seconds, and check whether 
P(2) = 1 in 1/4 seconds,. . . , and check whether P{n) holds in 2~" seconds, and 
so on. It is easy to conclude that if P find an answer, it will return the answer 
in one second, otherwise it will return false after a second. Considering the sum 
of geometric series, the proof is trivial. So the upperbound of the time for P to 
solve any question is 

Apparently, if P does exist, its power is extraordinarily great, for it can even 
solve Turing's Halting Problem in one second. 

So far we have seen two idea to implement the PLATO machine P. However, 
unfortunately, neither of them are successful. The first one is to construct the 
machine according to the definitions of H.Weyl. Apparently, it is difficult, for 
people do not believe that time is infinitely divisible. The second one is to make 
use of the theory of general relativity. However, the computing system will also 
exhaust the resource of the universe which make the observer cannot get the 
answer. 
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3.7.3 Recursive function whose derivative is not recursive 

April 1970, J.Myhill published his astonishing result[B]: There exists a recur- 
sive function, whose derivative is not recursive. In order to understand the 
principles of the construction, knowing the following fact about the recursive 
functions (whose domain is M) should be helpful. 

Theorem 3.15 Suppose f is a real-valued function, {/«} is a series of recur- 
sive functions, if there exists a recursive function e : N — >■ N s.t. Vx £ Ik > 
e(n) \fk{x) — f{x)\ < -^ , then f is recursively computable. 

J.Myhill's idea is to build a non-trivial structure(slope or bump) in the 
neighborhood of 2^" in interval [0,1], where n G £/, and jz/ is a recursively 
enumerable, nonrecursive set. Otherwise f{x) = 0. However, in order to make 
the function computable, the scale of the structure should shrink as the n is 
enumerated recursively, or rather, should be smaller than the bound in the the- 
orem above. As a result, the derivative of the function is intuitively hard to 
compute, and on the other hand we can proof that it is not recursive, because if 
we could compute it we can use the result to decide whether \x~\ is an element 
of £/ generally, contradicting the nonrecursiveness of .s^/. 

Specifically, suppose 

r a;(x2-l)2, if-l<x<l; 
^^ > -\ 0, if \x\ > 1. 

It is easy to verify that e{-l) = 9(0) ^ 6(1) ^ 09'{-l) = e'{l) = 06i'(0) = 
land 9min = 0{— 1/^/5) = —\9max = 9{+l\/b) = +A. We call 9 a bump of 
length 2 and height A. Then the function Oajiix) = {(3/\)9{x/a) satisfies the 
following conditions: 

Oc^i-a) ^ 9^pi0) = 9^pia) ^ 0, 9'^p{-a) = 9'^^{a) = 0, 9'^^{0) ^ 9/Xa, 

-I3<9ap{x) < 13 {a<x<a.) 

For each n e jz/we shall construct abump: 9a„i3„ at 2^" i.e. 
if n G ^, (5 e [-an, +a„], /(2-" + 6) = 6'„„"^„'((5), otherwise f{x) = 0. To 
make / well-defined, parameters a„, /3„, n S jz/ is defined as 

„ r, — k—2n—2 Q c\ — k—n—2 

where n = h{k) and h is a function enumerating £/ without repetitions (It is 
easy to proof that if there exists a recursive function enumerating jz/, then there 
exists such function with no repetitions). 

For physicists, does J.Myhill's results imply that if an object move under the 
condition that the displacement and the time satisfies the following relations 

rm = fm = / ^""'^"^'^)^" ^ '^^' if^ = 2"" + ^, Se[an,+an]; 
^ ' ■'^ ' \ 0, o.w. 

The speed v = r'{t) will be a physical quantity which is not computable? 
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3.7.4 Physical States which is not computable 

Pour-El et al published their results in 1997: for a differential equation, one can 
design a specific initial state to make the solution after t (t could be take some 
computable value)seconds is nowhere computablejS]. 
Consider the IVP of the following wave equation: 

{d^u ^ d^x , d^ j_ d^z 
u{x, y, z, 0) = /(x, y, z), |^(x, y, z, 0) = . 

where {x,y,z) £ M^, t e [ 0, +oo ) for aU / e "^^ this IVP has a form of 
solution known as Kirchhoff 's formula: 

M(^,i)= ff [f{~^ + tn)+t\If{~^ + tlt)-lt]da{lt) 

The conclusion Pour-El get is the following theorem: 

Theorem 3.16 For all compact set Z? C K'^ x [ 0, oo ), there exists a com- 
putable junction f{x,y,z) G ^^ , s.t. the corresponding solution u(x,t) is not 
computable in the neighborhood of any point in D. 

Pour-El et al construct the initial value through the uncomputable real number 

Apparently, one can conclude that in this wave equation, the initial state is 
computable but the state m(0, 0, 0, 1) is a state which can not be compute. 
For us, can we safely conclude that 

{Turing Computable} C {Physical Computable} 

but 



{Turing Computable} ^ {Physical Computable} 



? 



3.7.5 A few Comments 

In the above scenario, the use of (actual) infinity is their common theme. They 
ask the system to run for infinite steps or just encode the solutions into real 
numbers. It is easy to find out that adding either of these two assumes into a 
physical system will make the original system extraordinarily powerful. 

For example, we can throw a particle onto a plane [0, 1] x [0, 1] at ran- 
dom(obey the uniform distribution), then we can proof that with high probabil- 
ity, the x-coordinate(or j/-coordinate) of the center of the particle will indicate 
a non-recursive real number. In fact, in cell [0, 1] x [0, 1], the Lebesgue measure- 
ment for the recursive real numbers is 0, while the rest is 1, i.e. 

m([0, 1] X [0, 1] n Kr) ^ 0, m([0, 1] x [0, 1] n R^) = 1 

This is geometric probability and consider the uniform distribution, the proba- 
bility of the either event of the two are just their measurement. Therefore we 
can look the a;— coordinate as a function with respect to the digits. 

However, does the strict plane really exist in the physicla world? We just do 
not know. 
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We propose the some levels for / which is not computable. Suppose o is an 
operator, 6 is a physical implementation of o and operators A , V always exist. 
For f € F where all elements in F are non-recursive functions, we have the level 
of existence as follows. 

• Existence-I 

3p{A opoV = f) 

• Existence-I* 

3p3V3F{{f (eF)A (Pr{ a opoV ^ f\f e F} > 0)) 

• Existence-II 

3A3p(A o p o V = /) 

• Existence-II* 

3p3V{AopoV = /) 

• Existence-Ill 

3p3A3V(Aopo V = /) 

According to the levels we proposed above, assume the space is continuous, we 
can find out that P e Existence-II, J.Myhill's function / G Existence-I, Pour- 
El's construction (j> e Existence-II, our example x S Existence-I*. Apparently, 
we wish to get the examples in Existence-Ill. 

4 Physical Resource Complexity for Quantum 
Computation 

4.1 Physical Resource Complexity for Quantum Compu- 
tation 

For general quantum computation, we only need to explain the definition of the 
physical state set and the required evolution operators. More over, we only talk 
about Monte Carlo styled quantum algorithms. 

According to von Neumann's four postulates for quantum mechanics, we 
require that the state of any representation should be vectors in Hilbcrt space, i.e. 

fi CH 

and the evolution operators should be unitary, i.e. 

H e U{n) 

Without loss of generality, we can assume that the measurement operators 
is projection operators(POVM could be substituted by projection operators 
through adding more auxiliary qubits) 

4.1.1 The RCEF for Quantum Computation 

The resource cost by a computation is 

9l^(Ao-HoV) = (T(Ao-HoV), 
6(Ao-HoV), 
£(Ao-HoV), 
©(Ao-HoV)) 
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Our definition here is special a case of the one in Sec-III. Suppose our discussion 
is restricted to QCM, i.e. we have finite kinds of universal quantum operators, 
then the number of gates used and the depth of the whole circuit will be the 
main parameter which should be took into account. It is easy to find out that 
this definition is similar to that of quantum circuit model. One of the difference 
between them is that we will also take the cost of dcsign(usually this costs time) 
of a new circuit into account. Though in most cases, this will not cause great 
difference from the result given by QCM, however, we don't think we can safely 
ignore the potential exceptions just because it is usually easy to expand the 
scale of some circuits. 

So far, people always assume that qubit is relatively easy to prepared. At 
least in the asymptotic sense, no matter how difficult to prepared a quantum 
bit, the cost should be bounded by a constant. We will also do this. 



4.1.2 Deutsch-Josza Algorithm 

Deutsch-Josza algorithm is one of the most successful algorithms in the early 
years. The corresponding problem of the algorithm is: consider two sets of 
functions: 



A: ^1^ : {0, . . . , 2" - 1} ^ {0, 1}, Vx(^(x) = 0) 



B: ^1^ : {0, . . . , 2" - 1} ^ {0, 1}, {x\^{x) = 0} = {x\^{x) = 1} 



Apparently we have A n B = 0, now suppose / e A U S and there is an oracle 
to compute /. We are required to decide whether / G A or not. It is no doubt 
that people wish to reduce the times of query the oracle as much as possible. 

Note that the cost of implement the oracle is not taken into account, because 
we assume we have implemented it. 

The algorithm needs a trivial input -00 = |0)'^"|1), and used the gate H^"(E) 
H onto the state^'o a nd get V'l, i-e. 

^1 = (if»"®i7) (|0)^"|1)) 
75((|0) + |1))(0| + (|0) - |1>)(1|). By induction we have 



Note that H 



V2 



^^" = 4^E(-ir'i^)(2^i 



x,V 



where i.e. x • y = ^Xi Ajji. So we get: 



01 



H] (|0)«"|1)) 



iJ« 



= ^(^E(-iri^)(yi)io>^" 



|0)-|1) 
V2 



(by orthogonality) 



V2 



-h^EM 



x=0 



|0)-|1) 
V2 
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Now use the oracle Uf : \x, y) -^ \x,y ® f{x)) onto the state i/'i to get 11^2 



tp2 



/ 2"-l 



= t^'e'^X/w® 



x=0 




1 ''5^\_1)/W|^) 



|0)~|1> 
V2 



At last we use H®" ® / onto V'2to getip^: 



V'a = (h'^"®i){^ E i~^y^"''\^) 



x=0 



|o)-|i) 



V2 



= EE 



(-1)- 



'i-> 



2" 



|0)-|1> 



V^ 



The observer is supposed to check the first n qubits, note that the amplitude of 
|0)®" is E^(-1)^^"'V2"- If / e ^, fix) is constant and the amplitude of |0)^" 
is +1 or —1. So the amplitude of another cases should be zero and the observer 
wiU get 10)*^". On the other hand, if / G i?, the amplitude of jO)*^" wiU be zero. 
So the observer will always get a non-zero vector. 

In our opinion, the procedure could be written as follows. 



Hi 

^2 

A 



Initialize the state|0)®" (g) |1) 
Generate the whole circuit 

Uf 

JJfg,n ^ J 
l:\P^){P^\ 



Let •H = •Ha o ^^2 o -Hi 

(A o ■H o V) (|x)) = (a o -Hg o -^2 o -Hi o v) (\x)) = P{f e B) 

Though Deutsch-Jozsa Algorithm is great, someone still think it is not very 
useful. In addition to the fact that the problem they studied is not very impor- 
tant, there does exist an efficient classical probabilistic algorithm to solve the 
problem with high probability. 

4.1.3 Shor's Algorithm 

Shor's Algorithms for prime factorization and discrete logarithms pUl [TTj is so 
far the most exciting quantum algorithms. The appearance of Shor's Algorithms 
is the greatest challenge to strong Church Turing Thesis. 

Shor's Algorithms depends on a technique of so called "quantum Fourier 
Transform". But of course QFT is not enough. Shor's Algorithm is totally 
non-trivial and marvelous, and few people can produce any algorithms like that 
easily. 

In order to understand Shor's Algorithm, it may be enough to gain a clear 
idea of quantum ordering algorithm. This is the only subprogram in the Shor's 
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Algorithm which has to be implemented by quantum computers so far, and it 
is really the most important subprogram. 
First, note that 

r-l 

> exp(— 27risfc/r) = rSko 

s=0 

and define lu,) as follows 



r-l 
A V^ -2irisk /r 



fe'=0 



X modiV 



According to the fact above, we can get 



r-l 



-y 



^l-Ktskj r 



\Us) = 



s=0 



a:;''mod N 



In fact 



1 "'^ 



r-l 



g2^isfe/r| 



J_ Y^ / 2iTisk/r _}_ Y^ -2ixisk' jr 



s=0 



s=0 



fe'=0 
27risk /r 



x'^'modN) 



_. r — 1 / r — 1 

s=0 \ fe'=0 

1 v^ v^ f2Tris(k — k') 

-l.l.cxp' 

k'=Q s=f) 
1 '■"^ 



a;'''modiV\ 



x*^ modA^ 



k'=0 

= x'^mod A 
In particular, when k — Q, we have 



c'^'modAf ) 



r-l 



Us) - ID 



where L = [log(A^)] . 

Suppose Ux,N satisfies Ux.N\y) = |a;y(modA^)). Considering 1^*^^ and the 
fact that the permutation on orthnormal basis can be represented as a unitary 
operator, one can know for sure that Ux,n is unitary. What's more Us is a 
eigenvector oiUx,N , the corresponding eigenvalue is e~^ since 



1 '■"^ 

Ux,n\us) = -^y^i 



fc=0 



a;'=+VodA\ ^e^\us) 



Reverse the results above, we get the first half of the quantum ordering 
Algorithm, which complete the following task: 
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Il\»-^ ^ -— 



r-1 

El 

s=0 



,-l2" 



Modular exponentiation 




x'^'mod N 



where the state -h= 






•^'^**'^/''|us) is the one we desire. Apparently the eigen- 
value contains the information of r. So as to extract the information, we need a 
sub-progress named " quantum phase estimation" which based on inverse quan- 
tum fourier transformation. One can verify that if t is large enough, such like 
t = 2L + I + [log (2 -I- ^)] , for each s G {0, . . . , r - 1}, we will obtain the 
estimation of (p « s/r accurate to 2i + 1 bits with probability at least (1 — e)/r. 
Through the continued fractions algorithm, we will get r with high probabil- 
ity (According to PNT). 

In our opinion, the procedure above could be written as: 



V E 


E Initialize the State] 0)'»*(8) |1) 




Generate the whole circuit 


Ui = 


E iJ«5" (g) H 


H2 E 


= CU^,N 


^3 E 


E FTt ® I«'^ 


Ui E 


E CF(g) I'^^ 


A E 


- Y.\P^){PA 



It is easy to check that except the 7^3, all operators cost polynomial time with 
respect to log N . The complexity of operator modular exponentiation and con- 
tinued fraction are both 0{L^), which are two most time-consuming subproce- 
dure of the whole algorithm except the "Ha (inverse quantum fourier transform). 
Note that "Ha is indeed not an operator which could be implemented by 
polynomial universal gates. Consider a family of gates used in Tis which is 
usually noted by Rk{k G {2, . . . , L}) 



Rk 



1 







e2-V2'= 



In other words, the original Shor's Algorithm is not a algorithm with super- 
polynomial acceleration. In order to overcome this. Coppersmith created a new 
algorithm called the AFFT(Approximate Fast Fourier transform) jl2j which can 
substitute for the procedure QFT. 



4.1.4 Grover Algorithm 

Quantum Search Algorithm J13j. also known as Grover's Algorithm, is another 
quite successful quantum algorithm. Though this algorithm is not faster than 
the fastest classical search algorithms super-polynomially, one can proof it is 
the fastest one considering quantum mechanics. Therefore, the complexity of 
the algorithm is the complexity of the problem it deals with. 

The crucial subroutine of Grover's Algorithm is the Grover iteration, often 
denoted by G: 



Apply Oracle O : \x)\- 



(-l)/(-)|x)|- 
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• Apply Hadamard GatcsiiJ*^" 

• Perform a conditional phase shift(2|0) (0| — /) on the computer, with every 
non-zero bases receiving a phase shift of — 1 . 

• Perform Hadamard transformation H'^". 

Note that iJ^"(2|0)(0| - /)H«3" = 2|V')(i/'| - / One can proof that Grover 
iteration can be looked as a rotation in the plane spanned vectors which denoted 
the right answers and the wrong answers. 

Let E^ be the sum of all the vectors which indicate a solution to the search 
problem, E" the rest. Define normalized states: 

thus the initial state \tp) = ^^Ij^ '^x^ol^) could be represented as 



The action of Operator O is 0{a\a) + b\f3)) = a\a) — 6|/3), which could be 
looked as perform a reflection in a/3— plane. Similarly Operator 2 1 ■)/))( V'| — / 
also performs a reflection in a/3— plane. Thus both two reflections which could 
be looked as a rotation occur in the a/3— plane. Let cos 9/2 = ^/{N^^M)/N , 
s.t.\i)) = cos0/2|a) + sin0/2|/3), apply the iteration once makes |V') become 

G|-0) =cosy|a)+smy|^) 
k times use of Grover's Iteration will lead to the following result: 

G'^IV') = cos ( ^— ^ ) I") + sm ( ^~' 



Since 1^-) = yJ{N - M)/N\a) + y/M/N\l3), we just need to rotate \tp) 
arccos ^J MjN radians to the one which is parallel to vector |/3). So repeat- 
ing G for _R = [ J- ] times will get ji/i) to within an angle 0/2 < 7r/4 of 

1/3). This is a 'good' state, for people only have to repeat the experiment for 
expected constant times to get the solution to the problem (Consider geometric 
probability distribution: E[X] = 1/(1/2) = 2). 

Apparently R < [7r/20] , suppose M < N/2 then we have f > sin f = J^. 
Thus, we obtain: 

Ttt flV 

R< -\ — 

- 4 V M 



in other words we need repeat G for i? = 0{^N/M) times. 
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4.1.5 Quantum Simulation and Quantum Algorithm 

Quantum Lattice Celluar Automata(QLCA) and Quantum Gas Automata(QG- 
A) are two familiar ideal models in the research of quantum simulation [T51 . 
Meyer, Boghosian[T5l [TBI [17] have obtained their results respectively by using 
these models, that is, they construct some quantum algorithms which demon- 
strate exponentially speedup in such models. For Bohosian, the object they 
tried to simulated is a QGA which obey lattice Boltzman distribution, where 
arbitrary fields can be concerned with. They have proofed that the complex- 
ity of simulation is only related to the dimension of the lattice, but almost 
has nothing to do with the number of the particles. However, the number of 
particle always cause exponentially hardness on a classical computer. In fact, 
Boghosian's results imply that it is almost impossible for a classical computer to 
simulate one evolution step of a quantum system including dozens of particles. 

We've mentioned that it is the difficulty of quantum simulation that makes 
people believe quantum mechanics can provide enormous power of computation 
in the early years. 

Note that in this article we do not care about the hardness of simulations. 
Generally speaking, the hardness of simulation has nothing to do with the one 
of computation. For instance, people may find it difficult to simulate some 
classical celluar automaton according to the given regulations, however once the 
tedious work has been completed there often exists some more simple methods 
to produce the series. A typical example is that the regulations of an automata 
actually cause a circle with a finite period in the series. The same thing can 
happens to quantum simulations too. 

However, it is important to know that there must exists some cases in which 
simulations and computations are equivalent. These extreme cases often ap- 
pears when the length of regulations is near the Kolmogorov complexity(lower 
bound of description) of a series. Still, strictly speaking, at present no one 
can proof that polynomially universal unitary operators really cause exponen- 
tially difficulty in classical computation. To understand this, just consider an 
easy but helpful fact that almost all the problems we want to efficiently solved 
on a quantum computer are in the class BQP, and we have BQPCPSPACE. 
Unfortunately PSPACE=P is not totally impossible. Of course most people 
don't believe this is true, since this would imply that Shor's Algorithms can be 
polynomially simulated on a classical computers. 

Now we discuss how to extract a corresponding quantum algorithm from 
a method of quantum simulation, which is believed to be exponentially faster 
than any classical one of the same target. 

On a high level, we should do following things: 

• Find a family of experiments of quantum mechanics which can be effi- 
ciently simulated by quantum computers but are believed to be hardly to 
simulate and compute by classical computers 

• Design a 'good' problem about some non-trivial properties of the last state 
of the system, which makes quantum computers able to present the answer 
to the observer quickly. 

Designing the problem is a crucial step. In most cases, though we may have 
quickly obtained the probabilistic distribution very close to the real experiments. 
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we can not know the whole information in short time. So first we have to ask 
a question which can be easily verified by any quantum computers containing 
the whole quantum information of the system. 
For example, we can ask a question such like: 

• What the number of the tiq— th digit of the probability of a certain system 
arriving in il'(i7' C 17)? 

The problem of this method is that in high dimensional spaces, it is very 
likely that the probability of the set i7' is exponentially close to zero, which 
actually enables a classical computers to guess zero without running and get 
the right answer in most cases. 

Now we propose our version: Suppose is the wave function of the system 
we've simulated and \(p{X)\'^,X C il is the probability of x appear in X. Try 
to find two subsets A, B C fl s.t. 



\^{AW 



3 

7- |</-(S)P 



< 1 



and determine the value of the no— th digit of (f){A). 

For the systems which (probabilistic) Turing Machine cannot simulate in 
polynomial time, the question above is intuitively hard to answer, though up 
till now no one can proof or disproof it. 

On the other hand, if these systems can be efficiently simulated by quan- 
tum computers, repeating following procedure will ensure us to find the answer 
relatively much faster than any probabilistic Turing Machine of the same aim. 

Definition 4.1 (Vector of normal vectors x) 





Vi/ 



1 



, . . . , 





voy 



Definition 4.2 (Procedure Pq) Pqiln pseudo-code): 
while(find the answer) 

{ 

Mid-cut the space il by super-plane whose normal vector is Xi 

Suppose the two spaces is i7i and 0,2 

if(the condition is satisfied(verifies by testing)) 

{ 

halting 

} 
else 

{ 

fl ^ min|0|{r2i,fi2} 
i + + 
} 
} 
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5 Conclusions and Future Works 

We formally proposed the theory of physical computation, define the concepts 
of resource and complexity. Several examples, including classic mechanics and 
quantum mechanics, were discussed and analyzed under the framework of phys- 
ical computation. A technique, which is used to converse a method of quantum 
simulation into a quantum algorithm, is discussed. 

This is a exciting field, we believe there is more exciting topic to discussed. 
A very interesting question is: can we find a physical mechanism as the fastest 
implementation of an arbitrary functions? 

In Sec-III, we talk about the question of calculating the centroid of an ob- 
ject. We thought it is the limitation of dimensions (only three dimensions) hide 
the advance of the method we mentioned. We conjecture that this method 
has a excellent counterpart in high dimensional cases. We'll have a try in the 
(quantum)statistics mechanism. 

In Sec-IV, we talked about quantum simulations and how to construct a 
clever problem to induced a quantum algorithm. Actually, we conjecture that 
the problem we construct is a hard one in class #P, for these questions have a 
counting style. However, we are not sure about whether the designed questions 
could be in #P — hard under some specific statistical models. We shall try to 
work on this in the future. 

We've mentioned that we assume that polynomial qubits is polynomially 
hard to prepare. However, it is harder to control the qubits as the number of 
them increase[TSj so far. So one can still conjecture that preparing qubits itself 
is a " complicated computing" , and the results up till now can be explained as 
someone displace the resource consuming procedure, just like DNA Algorithms. 
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"<) . Inspired by the work of Feynman, Deutsch, We formally propose the 

theory of physical computation and accordingly, the physical complexity 
' (-H \ theory. To achieve this, a framework that could be used to evaluate al- 

O ,1 most all forms of computation making use of various physical mechanisms 

is established. Here, we focus on applying this to Quantum Computation. 
As a preliminary study on more general problems, some examples of other 
physical mechanisms are also discussed in this paper. 
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1 Introduction 

1.1 Quantum Computation 

The research of quantum computation has been lasting for about 30 years since 
R.Feynman proposed the concept of so-called 'quantum computer' in 1982 2]. 
Founding out that there do exist some quantum systems which are suspected 
cannot be efficiently simulated by classical computers, early researchers natu- 
rally speculated that quantum mechanism itself may provide stunning power 
of computation. In order to strictly define what is 'quantum computation', 
researchers introduced various new computational models, including Quantum 
Turing Machine [1] and Quantum Circuit Model[5]. However, at that time, 
no convincing evidence was discovered to support the conjecture that quan- 
tum mechanism can really be used to speed up the computation of some hard 
problems greatly. 

D.Deutsch found the first evidence that quantum computers may surpass 
the Turing Machine £ in query. In fact, he constructed a special scene in which 
DTM has to query the oracle for 0(2") times to find the correct answer(for 
certain) in worst cases while QTM need just once query in all cases. 

One of the most remarkable results is quantum factorization, which is due to 
Peter Shor pUUTT] . The best classical algorithm for factorization so far has to run 
for 0(log L) steps. However, Shor showed that one can use a family of quantum 



circuit, which contain 0(log^ L) gates and needs only 0{L^ log L log log L) op- 
erations to get the right answer. 

Grover's Algorithm ITS] is another successful example of quantum algo- 
rithms. This algorithm can be used to search a database without structure. 
It is easy to proof that the time complexity of this problem for Turing Machine 
is 0{n). However, there does exist a quantum algorithm whose time complex- 
ity is 0{y/n). Since it has been proved that this is the optimal algorithm for 
all algorithms that considering quantum mechanics |18j , so the complexity of 
Grover's algorithm can be looked as the quantum complexity of this problem. 

One of the most important reason that why Quantum algorithms(especially 
Shor's algorithm) seem so interesting to many computer scientists is that their 
existence indicate a huge challenge to strong Church- Turing thesis, which states 
that Any model of computation can be simulated on a probabilistic Turing ma- 
chine with at most a polynomial increase in the number of elementary operations 
required. 

1.2 Physical Computation 

On the other hand, with the exciting research in quantum computation as 
well as other new paradigms of computations(e.g.DNA computation), the idea 
that we may just look physical processes as computations(not just the Tur- 
ing Machine)was also developed. The seeds of this idea can be traced back to 
Feynman[2 , Deuthsch[3, and Pitowsky[7j et al. 

It is not very hard to understand and appreciate this idea, for at first glance, 
this point of view has at least three benefits: 

• It can include the concept of classical algorithms easily, for an algorithm 
on Turing Machine(its physical implementation)can also be thought as a 
family of physical processes (and the corresponding measurement). 

• We can try to solve some special kinds of problems with less time or space 
than the lower bound with respect to Turing machines. 

• Being the ones which could be directly simulated, some physical methods 
can also enlighten us to design smart algorithms on Turing machines. 

What's more, currently, it seems that we cannot exclude the possibility that 
there does exist a family of physical processes which can help us to calculate 
some problems which cannot be solved by a universal Turing Machine in prin- 
ciple. 

However, because of vagueness and extraordinary generality, the theory of 
so-called 'physical computation' has a significant defeat yet. 

• In many cases, people cannot decide how to define the resource for a 
'physical algorithm'. And as a result they cannot proof or even conjecture 
formally whether a 'physical algorithm' is really superior to any algorithms 
of TMs with the same cxtensionality. 

Note that the theory of quantum computation is almost free from such defeat, 
for researchers have completed the formal definition of the computational model 
of quantum computation in the early years. Roughly speaking, things tend to 
go wrong when: 



• People adopt a design on which the physical postulates it depends is just 
an empirical one. 

• More than one different systems of physical postulates are used. 

1.3 The structure of the article 

The structure of this paper is as follows. In Sec. II we shall introduce two 
well-known models of quantum computation and the definition of complexity 
respectively. And after that we will formally establish the theoretical foundation 
of physical computation and propose the theory of physical computability in Sec. 
III. In the beginning of Sec. IVwe try to use the theory of physical computation 
to reanalyze the quantum algorithms. In the end of Sec. IV, we focus on 
the topic about how to construct problems which take advantage of quantum 
simulations. 



2 Models of Quantum Computation 

2.1 Quantum Turing Machine 

Quatnum Turing Machine was first introduced by BeniofF[T] in 1980 and was 
developed by Deustch and Yao. The modern definitioin was given by Bernstein 
and Vazirani in 1997 4 . 

Definition 2.1 (Quantum Turing Machine, Bernstein 1997) Let C be a set 
of complex nmber a satisfying: For each a, there exists a polynomial time 
algorithm to compute the value of Im{a) and Re{a) close to 2~" within the 
true value. 

A Quantum Turing Machine M is defined as the triple (E, Q, 6), where S is 
a fintie alphabet with an identified symbol ^ , Q is a finite set of states with 
an identified initial state qq and final state qf ^ qo; 6, the quantum transform 
function 6 : Q x Y. ^ C^'"^''^^'^^ The QTM has a two-way infinite tape of 
cells indexed by Z, and a single read/ write tape head that moves along the 
tape. We define configurations initial configurations and final configurations 
exactly as for DTMs. Let S be the inner-product space of finite complex linear 
combinations of configurations of M with the Euclidian norm. We call each 
element phi € S" a superposition of M. The QTM M defines a linear operator 
Um : S ^ S called the time evolution operator of M as follows: If M starts in 
configurations c with current state p and scanned symbol a. The after one step 
M will be in superposition of configurations ip = ^^ aiCi, where each non-zero 
ai corresponds to a S{p,a,T,q,d), and Ci is the new configuration that results 
from applying this transiton to c. Extending thi map to the entire space S 
through linearity gives the linear time evolution operator Um- 

Definition 2.2 If Um can keep Euclidian norm, then we say M is well de- 
formed. 

Theorem 2.3 If QTM is in the superposition tp = '^- aiCi and is observed, the 
probability of the observer gets the configuration Cj is |aip, and then M is in 
the state -0' — ci. 



Theorem 2.4 A QTM is well-deformed if and only if its time evolution oper- 
ator is unitary. 

In QTM, the number of the read/write tape head moves during a computation 
is the cost of time. 

Theorem 2.5 There exists a universal QTM, which is polynomially equivalent 
to any QTMs. 

2.2 Quantum Circuit Model 

The first quantum circuit model was due to Deutsch. Then quantum circuit 
model was improved by Yao[5], who also proved that for any QTM, there exists 
a uniform family of quantum circuit which is polynomially equivalent to that 
QTM. 

Not like QTM, quantum circuit model tends to describe an algorithm by us- 
ing universal quantum gates and circuits without loops. Quantum circuit model 
does not need infinite many quantum gates, but finite many quantum gates 
which called the universal quantum gates. It has been proved that Hadamard 
Gate, phase gate, C-NOT Gate and tt/S Gate are universal. For any finite 
dimensional U operators, we can always approach it effectively by means of a 
universal family of circuits ^ , which only consists 4 gates above, i.e. 

Ve(3ne '^),^([/,C/„) = max||(C/-J7„)|V')|l <e 

W 

The scale of a quantum circuit is defined as the number of the universal 
gates and the depth is defined as the longest path from input to output, if the 
gates is looked as a vertex. 

Both Quantum circuit model and QTM are important models of quantum 
computation. But we do not know whether they are the most natural models of 
quantum computation or do they fully take the advance of quantum mechanics, 
no matter in the theory of quantum computability and quantum computational 
complexity. 

3 The Theory of Physical Computation 

3.1 Observer 

Measurement is in terms of observer. Though there are many differences among 
people's opinions about the exact definition of human beings, we prudently 
assume that an observer is classical, that is, the observer will never get incom- 
patible results during one measurement. 

In this article, we will never use the terminology such like 'a observer of the 
observer', or in other words, by 'observer' we always mean the last one outside 
the whole experiment. 

In order to unify various forms of results, we require that the observer only 
accept the symbols on a tapeQust something like the one of TM) and also only 
use this to initialize an experiment. 

So we define the legal inputs and outputs as the elements in set E+, where 



S = {0,1,*,.} 

and S+ the finite string composed by elements in S. 
The concept of observer is crucial to our theory. 

3.2 Physical States 

We use (usually finite) attributes which may contribute to the computations to 
label the physical states. In addition, though may not be actually concerned 
in every computation, three fundamental quantities, namely, space, energy and 
mass are always included in a state for the sake of analysis of resource and 
complexity. 
So we have: 

n C {xi}^' X {xa}-^^ X • • • X {xn}-^" X {m}®' x {s}® x {e}*^ 

Or more gcnerally(Quantum) 

n C {Xil^i X {X2}'^' X • • • X {Xn}"^- X {C"}®' X {C}® X {C}'^ 

For simplicity, fundamental attributes are usually omitted, i.e. 

n C {Xi}^^ X {X2}^' X • • • X {x„}^" 

For a certain attribute Ai , what really matters is its type which is constrained 
by its dimension. Note that dimensionless quantity(e.g. friction coefficient) can 
also be assigned to a null type. When a quantity is expressed by other quantities' 
combination, it's dimension type should be preserved, or rather, any equations 
should be dimensional balanced. 

For example: 

^ID:ML''T-''] .._ ^[D:M]g[D:LT-^]f^[D:L] ^ ;^g/j[^:AfL^T-^] 

are all dimensional balanced. 



3.3 Physical processes and the operation 



o 



Physical process on a state space il is a set of physical state whose elements are 
labeled by moment t{t e [0,T],r e R+). 

P = { (T,P)| reM+,P: [0,T] ^ 17 } 
If two physical processes on Vt satisfies 

{t12Pi){t^iPi) - (^2P2)(0) 

we can define operator o : {P} x {P} -^ {P} i.e. 

^2 o Pi - Pa 
satisfies: 



1. 7riP3=7riPi+7riP2 

2. if < t < TTiPi (ttsPs) W - (^2Pi)(t) 

3. if TTlPl < t < TTlPl + TTlPz (^2P3) W = (^2P2)(t - ^iPl) 

For convenience, we introduce 1>P<I as the initial state of P, and <Pl> the final 
state of P, i.e. 

C>P< = (^2P)(0), <P^ = (7r2P)(7riP) 

3.4 Physical Operator and the operation of operator 

Physical operator is a tuple whose first component is a state x in fl and the 
second component is a physical process whose initial state is x, i.e. 

= {{x,P)\xe Q,\>P<i^x} 

In particular, a deterministic physical operator O means:0 is a physical opera- 
tor, and 

if TT20{xi) ^ 7r2 0(2:2), then xi ^ X2 

if we only care about the effect the operator do to the initial state, we can look 
operator as a mapping in f2,i.e.O : fi — >■ O. 

The operation between two deterministic physical operator is defined as 
follows (suppose Oi,02 are productive): 

03 = 020 Oi 

which satisfies 

yxiO^ix) = O2(<]0i(a:)o) o Oi(a;)) 

In some more general cases, it is useful to talk about non-deterministic phys- 
ical operators or random physical operators. A random physical operator O 
contains the tuples which has the same initial states but different physical pro- 
cesses, i.e. 

d = {{x,P)\xe Q,>P<i^x}. 

People cannot decide the output 0{x) just by the initial state a: e f7. 

Similarly, if just care about extensionality, we can look operator as a rela- 
tionship on il i.e.O : Jl x fi 

We can also define operations between two non-deterministic operators, if 
some preconditions are satisfied. To do so, we first expand the definition of 
some symbols. 

0{x) = {P\[>P<i=x} 

0{X) = {P\t>P<ieX cQ} 

<lO(x)[> = {y\y e rj,3P e 0{x)s.t. <] P[> = y} 

So O2 o Oi (if they are productive) can be defined as 

O2 o Oi(x) = {P2 o Pi|Pi e Oi(a:), P2 e 02(<Oi(x)[>)} 

Note that of all the processes created by random physical operators, their 
last states shall be exposed to outside world in the end by default. 



3.5 Physical Computability 

3.5.1 Deterministic Physical Computation 

Definition 3.1 (Deterministic Physical System) Deterministic Physical System 
^ is a Five- Tuple 

^ = (r!,I],V,H,A) 



where: 



• 



• 



S = {0,l,*,.}S+is the collection of finite string formed by elements in S , 
rix;+ is the the set of physical implementation of E+. 

il = {ipi,i E A} n ^ fl^+A is an index set il is a set of distinguishable 
physical states(labeled by their attributes). 



• V : ^^+ -^ reinitialization operator 

• H : fl ^ fl evolution operator 

• A : fi ^- J7 5]+ Measurement operator 

Since Hilbert's 6th problem has not been solved yet, i.e. the whole theory of 
physics has not been axiomatized, we do not know that whether there exist 
some additional fundamental mathematical constraints should be included in 
this theory, though Beggs et al. have discussed the axioms of measurement 
based on Hempel's axioms [T5]. Now, maybe the only restrictions here are the 
finiteness of the resource cost by a physical process and the finiteness of the 
attributes used to label a physical state set. 

As a result this system may looks looser than many classical computational 
models and may contains the ability to surpass all these models. We would like 
to let physicists to add more necessary restrictions into the system. 

Of course, when it comes to a specific branch of the physics, we can always 
know what is a legal states and processes. However, we wish to keep some 
freedom, i.e. to let the observer combine various axioms in physics so as to 
optimize the computations. 

Definition 3.2 (Partial Physical Computable Arithmetic Functions)For any 
partial arithmetic function / : N — > N is said to be partial physical computableii 
and only if there exists a Deterministic Physical System 

which satisfies 

If a; G don2(f), then, 

(Ao-HoV)(a;) = /(x) 

Similarly, we can define Total Physically Computable Arithmetic Functions 

Definition 3.3 (Total Physically Computable Arithmetic Functions) For any 
total arithmetic function / : N — > N is said to be Total Physically Computable, 
if and only if there exists a Deterministic Physical System (fi, E, V, "H, A) which 
satisfies 



Vx G Nwc have: 

(Ao-HoV)(.t) = /(x) 

In order to extend the definition of physical computability to non-arithmetic 
functions, we should take into consideration the precision of the measurement 
and computation. Therefore, we need a distance function to measure the pre- 
cision of two values and define the computability as the ability of computing in 
any desired precision. 

Definition 3.4 (Partial Physically Computable Functions) Given a partial func- 
tion / : A ^ Band a metric T) : B x B ^ M/ is said to be partial physically 
computable with respect to the metric T) , if and only if for any e > there 
exists (il, S, V, "H, A)s.t. for any x € Awe have 

if X G dom{f), 

p((Ao-HoV)(a;),/(x)) <e 

Similarly, we can also define Total Physical Computable Functions. 

Definition 3.5 (Total Physically Computable Functions) Given a total func- 
tion / : A ^ Band a metric T) : B x B ^ M/ is said to be partial physically 
computable with respect to the metric 2? if and only if for any e > there 

exists (rj, S, V, "H, A), such that, 

MxeA 

p((AoHoV)(x),/(x)) <e 

3.5.2 Non-deterministic Physical Computation 

On the other hand, many physical processes are considered to be non-deterministic, 
which enable us to implement so-called 'randomized algorithms' and 'quantum 
algorithms'. Our Probabilistic Physical System is defined as follows. 

Definition 3.6 (Probabilistic Physical System) Probabilistic Physical System £P* is 
a five-tuple: 

. where, 

. S = {0,1,*,.}. 

• n^{i:^,ie A}. 

• V : il^+ -^ riwhich is also called initialization operator 

• H* : ilxft which also called evolution operator, which is non-deterministic. 

• A : il -^ nY:+ , which is also called measurement operator. 

Non-deterministic does not necessarily cause probability, but let's convention 
that in this article we always discussed the randomness which has a probabilistic 
distribution. 

Definition of the computable functions by means of ■^* is an analog to that 
of ^. As an example, we define Total Non-deterministic Physical Computable 
Functions. 



Definition 3.7 (Total Non-deterministic Piiysical Computable Func- 
tions(Las Vegas)) For any total function / : N ^ N is said to be total 
non-dctcrministic physical computable function, if and only if there exists a 
five-tuple 

(f7,E,V,H*,A) 

s.t. 



Vx((Ao-H*oV)(a;) = f{x)) 

Because of randomness, for any identical inputs x, the system may call 
different process to compute. The above definition is the counterpart of the 
definition of the so called Las Vegas algorithm. 

Definition 3.8 (Total Non-deterministic Pliysical Computable Func- 
tions(Monte Carlo))For any total function / : N ^ N is said to be total 
non-deterministic physical computable, if and only if there exists 



(r!,E,V,H*,A) 



s.t. 



VxeN 

Pr{i<^ueni(x)occurs} > 2/3 

where 

Event{x) = ({A oH* o V){x) ^ f{x)\ 

3.6 Estimation of the Complexity of Physical Resource 

For the physical systems defined above, we can even ignore that whether there 
exists a physical mechanism in reality to implement it. Any functions which 
could be written as the composition of the three operators would be considered 
as computable (deterministic version). 
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But any experiments which implement a certain system will cost resource. 
We will focus on four kinds of resource, namely, time, space, energy and mass. 

One can also include 7 fundamental sorts of attributes considering 7 dimen- 
sions in SI(Systeme international d'unites). However, actually our decision is 
not a totally empirical one. In fact, many differential equations which are es- 
tablished to describe various phenomena involving these attributes are always 
related to energy, time, and space. After properly choosing unit to make all 
the constant into 1, we can just rewrite these dimensions in the forms of the 
combinations of 4 fundamental dimensions, which enable us to continue to use 
the 4 attributes to represent the increase of the resource. The counterexample 
may occur only when m or more than m new attributes appear in n equations 
and m > n, however these cases tend to be unlikely to happen, if these equa- 
tions(theory) we discussed are assumed to be 'enough' for some phenomena in 
nature. 



10 



Definition 3.9 ( Metric for Resource ) The resource of a physical process $H 
includes: 

• T: The (expectation of the)total time the whole process consumed; 

• &: The maximum of (the expectation of)the space the whole process con- 
sumed; 

• QJt: The maximum of (the expectation of)the mass the whole process con- 
sumed; 

• 2;: The maximum of (the expectation of)the energy the whole process 
consumed. 

and$R= (T,6,OT,€) 

In the above definitions, the metric of them could be selected as the common 
ones. Today, most physicists tends to believe that mass and energy are not 
independent, neither do time and space. But for convenience, we still focus the 
primitive forms of resource, for actually we don't care about the independence 
here. 

In the above definitions, we don't talk about the potential possibility that 
even time could be reused. 

We convention that the resource is with respect to an inertial system, i.e. 
the observers obtain their results when they are in an inertial system to the 
system running the 'algorithms', so as to rule out the paradoxes because of the 
theory of relativity. 

In many cases, we just cannot get a infinite precise estimation about the 
resource, but for our purpose, we actually do not need such things. Of course, 
there may exist some cases when we could not get an estimation without any 
promise of any precision, however, we will not use such processes to construct 
our implementation. 

Suppose the projections of the fundamental attributes (resource) are TTtxniS), 
7rg(S') and 7rg(S') Then the resource a physical process consumed is: 

TP = TTlP 

MP = max{7rOT(S'),5e Ran(7r2P)} 
&P = max{7re(S'),5eRan(7r2P)} 
£P = max{7re(5),S'e Ran(7r2P)} 

In general cases, when we have to discuss the process of superposition, the 
resource can be defined as: 

IP = TTlP 

MP = max{E[7rOT(S')],S'eRan(7r2P)} 
eP = max{E[7re(S')],S'eRan(7r2P)} 
€P = max{E[7re(S')],S'eRan(7r2P)} 



So it is easy to see that 



TP20P1 ^ 7riPi+7riP2 

mP2 o Pi = max{ mPi , srnP2 } 

&P20P1 ^ max{ mPi , 97IP2 } 
(£P2 o Pi = max{ mPi , 97IP2 } 
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According to the definition above, the resource of a non-deterministic physical 
operator O which is initiaUzed by a; € fi should be defined as: 

TO(x) = E[%Oi{x)],Oi{x) eO{x) 

mo{x) = E[mo,{x)],Oiix)eo{x) 

60{x) = E[eO^{x)],0^{x) eO{x) 
<BO{x) = E[e0^ix)],0^{x)eO{x) 

So for operators' operation, we have: 

%02oOi{x) = TOi(a;)+T02(<lOi(a;)l>) 
m02oOi{x) = max{9JlOi(x),OT02(<Oi(x)[>)} 

eo2oOiix) = m&x{eOi{x),eo2i<Oiix)\>)} 

(B02oOi{x) = m8ix{eOi{x),e02{<\Oi{x)c>)} 

3.6.1 Framework for the Complexity with respect to General Phys- 
ical Resource 

Definition 3.10 (Resource(deterniinistic))A resource the physical process 
which complete the whole computation consumed 91^ including: 

T<?. = T((Ao-HoV)(a;)) 
6^ = 6((Ao-HoV)(a;)) 

m^^ = m{{Aonov){x)) 

e^ = (£((Ao-HoV)(a;)) 

i.e. $H,3a = (1,32, S, 32,971^, €,32) 

Definition 3.11 (Resource(Las Vegas))A resource the physical process which 
complete the whole computation consumed d\^ including: 

T<32 = T((Ao-HoV)(x)) 

652 = 6((Ao-HoV)(a;)) 

9K52 = mi{AonoV){x)) 

^S^ = (£((Ao-HoV)(a;)) 

i.e. *H,^ = (1,32, 6^,9Jt,32, c£^) 

Definition 3.12 (Resource(Monte Carlo))A resource the physical process 
which complete the whole computation consumed $Hcs2 including: 

T<32 = T((Ao-HoV)(a;)) 
6^ EE 6((Ao-HoV)(a;)) 

ms^ = miiAonov){x)) 

(£<5» = (£((Ao-HoV)(a;)) 
i.e. D\gi = (T,32, (3,32,2H,32, C;^) 
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The corresponding concept of complexity should be defined as the resource 
consumed with respect to the length of the input. 

Definition 3.13 Complexity The complexity of a kind of resource is a func- 
tion of the length of the input x n = [log x~\ 

Conplexity^ (n) — max{T^(j,)|n — 1 < logo; < n} 

Conplexity„(n) = max{9Jl^(2,')|rt — 1 < logx < n} 

Conplexity^ (n) ~ max{(3^(j:)|n — 1 < logx < n} 

Conplexityg (n) = max{€^(a;)|7i — 1 < logx < n} 

Note: actually, the finiteness of resource cost is a necessary precondition for 
all "uniform" physical computation models. Another necessary precondition is 
that all the physical attributes should be at some trivial states (very easily to be 
constructed, e.g. 0°C, Om/s) before the experiments. If these requirement are 
not satisfied, the model will be a non-uniform one. This case is also discussed 
in details by Beggs et a/./7 ^ .FS ^ 

3.6.2 Some Common Examples 

It is interesting to find some new methods to compute problems without the 
help of universal Turing Machines. Through the ages, people have found a lot 
of such examples, the most famous of them are: 

• Measure the volume of an object by putting it into the water; 

• Obtain the centroid of an object by two suspension method; 

• Compute function sine by analog circuit; 

• Decide the path of minimum cost using Fermat's Principle; 

• Calculate the mean of numbers by the second Law of Thermodynamics. 
Actually, we can give even more similar examples: 

• By making use of resonance, we can easily find the desired tuning fork 
from a heap of them. Otherwise, we have to look up the label of them one 
by one and even have to compute the frequency one by one if there is no 
labels on them. 

• We can compute the square root of an given number x by the law of free 
fall. Prepare a vacuum tube T of length x and let it stand vertically, then 
let an object o which is small enough fall. Get the time t when it touch 
the bottom, and we have \/x = t/c, where c = (2/(7)^/^. 

• We can sort a scries of numbers through dangling poises by strings, where 
the strings satisfies Hooke's law. Given an array of numbers {xi} construct 
or find poises whose mass is just Xi , then dangle them by strings with the 
same stiffness coefficient. When the system is stable, the position of the 
poises with respect to their weight just indicate the relationship desired. 

However, it is hard for us to estimate the cost of the methods above just 
after we describe them informally. So we select a part of them to analyze next. 

Conventions: x is the representation of number in digits, [x] is the value of 
a;, \x\ means the attribute A has the value x. [x]^ means the representation 
of quantity x though not on the tape. 
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3.6.2.1 Mean of Three Numbers 

Given three numbers, compute the mean of them making use of law of thermo- 
dynamics. This idea comes from Pitowsky [7]. 

The strict description of the problem: Given three numbers xi, X2, X3 G 
[ 0, 100 ], compute 

X = (Accurate to two decimal places). 

o 

Pitowsky suggests that since all of the three numbers less than 100 and bigger 
than zero, note that the freezing point of water is C° and the boiling point 
of water is 100 C° under the one standard air pressure. So for each number Xi, 
we can prepare the corresponding water of volume V and temperature of XiC° . 
And then pour the water of three vessels into a bigger one, whose volume is 
V'{V' > 3V), and wait. After the water arrived at the balance point, measure 
the temperature. Of course, we assume that during the whole procedure, no 
calory is lose. 

Apparently, the physical state the method above deal with is the temperature 
of water, so we have 

n = {t\ti G [0, 100], i = 1,2,3}^. 

on the other hand, wc suppose the water is heat up from C°, i.e. the initial 

state of the experiment is ([0]-^, [0]-^, [0]'^). 

Therefore, the process could be depicted as following: 

V : S"*" — >■ O, heat up the water to the desired temperature 

V(a;i,a:2,a;3) = {[xi]'^ , [0:2]^, [xa]^) 

"H : rj — >■ fi, admixture the water of different temperature, the second law of 
thermodynamics is used 

?^([xiF,[x2F,[x3F) = ([sF,[sf,[xn 

A : f7 — > S"*" , measure the temperature of the water 

Ai[xf,[xf,[xf)=x 

For this problem, since the precision is finite, and there are only constant 
(three) numbers and the numbers are bounded, we can easily deduct that D\,^ 
is a constant. As a matter of fact, for Turing Machine, we can also find a 
constant resource costing algorithm which is just looking up a finite list to solve 
the problem. 

3.6.2.2 Sorting Without Repeat 

Description of the Problem: 

Input: Finite number series of length n: 

A = {x^\xi G Z+ n [0, Af](0 < i < n)}- 
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Output: Finite number series of length m: 

B = {xj\xj e ^(0 <j< m)}, 
s.t. if ji < J2 then Xj-^^ < Xj^ . 

Our plan is: for the given series, select a series of poises of length n, s.t. the 
mass of the ith poise is equivalent to the ith number. Dangling the poises from 
right to left by strings, whose restoring coefficient are k. Wait until the system 
is stationary, open the parallel light source and measure the projection onto the 
vertical ruler at the right end. The measurement could be done by machines 
and present the results onto the tape for observer. For Example, we can embed 
some photoconductivc diodes in the ruler by graduations, diodes who is not 
triggered should be read. 

The physical state the method is primarily concerned with is the mass of 
poise M, the horizontal positions of the poises X and the vertical ones Y, the 
projections Y' and the boole value B indicating which diodes is triggered, i.e. 

So we have 

V : S+ — > r2 (Select poises) 



v{(b7x,) = {®^[x,]''\2fm'' ef [(j,o)]^ 

"H : ri ^ ri (Dangle poises) 



n 



:x 






J 



]X 



[[x,]g/k]^ 

V®f[0-,o)] 



Y'xB 



W -M^ f7(0pen the parallel light) 



n' 



( mx^r 

[[Xr]9lkV 
V®,[(j,0)]^'x^ 






M 



[[x,]g/k]^ 
V®f[(j,eA(j))] 



Y'xi 



J 



A : ri — )■ S+(read the projection) 



/ 



M 



®r(FiJ 
[[x,]g/k]^ 



(xf) 



:jji yxji 



V®f[(j,eA(j)r ^^) / 



satisfies if ji < J2 then 

Considering the ideal implementation, we conclude that the Uig^ is linear, 
which is superior to Turing Machines using comparisons, for the complexity for 
them is proofed to be O(nlogn). However, there does exist Turing Machine, 
which is not based on comparisons, also has a linear time cost. 

Note that if the number series is boundless, the complexity of the method 
above will be exponential. This is the common defeat of most analog computers. 
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3.6.2.3 Volume of irregular shape 

For this issue, we shall restrict the range of the saying 'irregular' so as to rule out 
the objects with infinite length of description. So actually, we tend to discuss a 
subset of the set of all cases. 
Description of the problem: 

Inputs: point scries of length n:{xi,yi)(l < i < n), satisfies c + r < Xi < 
a — c — r, c + r < yi < b — c — r 

Outputs: The volume of the box of length a and width b and height ho, not 
including the series of cylinders(radius: r height: ho) which arc induced by 
the series of points. 

Our plan is simple. Assume we have a box of material of dense p, and a punch 
to extract circles from it. Then we measure the mass of the rest then divide it 
by its dense or just put it into water. V : S+ ^ il 

V(®r=ia(^^,2/^)) - [pho{A-m{uUc,))r[0f 

= [pho{A - m(ur^iC,))]^[p/io(A - m(ur^ic,))]^' 
-H2 : O ^ 17 

H2[ph„{A - m{U^^,c,))]^[pho{A - m{U^^,c,))f 

= [phoiA - ™(ur=iQ))]^[/.o(A - m(ur^iQ))]^' 

A : O ^ S+ 

A[pho{A - m{U^^,c,))f'[h„{A - m{U^^,c,))f ^ ho{A - m{U^^,c,)) 

Apparently the resource complexity for this method is linear with respect to 
the number of the points. However, because most people think that we cannot 
do infinitely measurement during one experiment, this method can only provide 
the result of finite precision. This is a good news to Turing Machines because 
this implies there exists a Turing Machine which is almost equivalently efficient. 

This may be astonish to someone, who may thought that a TM should at 
least solve the equations first. However, because of the finite precision, Turing 
Machine can just split the object into lattice and use the so-called scan-line 
algorithm to find the answer. 

3.6.2.4 The centroid of Irregular Shape 

Just as the last example, we restrict our topic into the same subsets of all cases. 
Description of Problem: 

Inputs: point series of length n:{xi,yi){l < i < n), satisfies c + r < Xi < 
a — c — r,c + r<yi<b — c — r 

Outputs: The centroid of the box of length a and width b and height ho, not 
including the series of cylinders which is induced by the series of points. 
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The method we suggest is similar to the last one, the difference of them is 
that this time we will record some points. 

V : E+ ^ rj 

Hi : 17 ^ fi (Suspend the box by Vq) 

= [phoiA m(ur.ic,))]*^[f— ^ + Vof[0f[0f" 

\c-Vo\ 



T-L2 : J7 ^ O (Suspend the box by V^) 



Mr ^ ^1 TAiSrnlS'rnlS" 



H3 : ft -^ ri(Extend the unit vectors: Get the point of intersection) 
H,[pho{A - m{iJUc^))r[-|^^^ + ^ol^^Iy^-T^ + V^f[of' 



C-^O , ,,iErC-Vo' 



Mi '^ "^1 I T/iEr'' ^0 , T/nS'rjS' 



c ■ 



A : fi ^ S+ 

AlphoiA - m(ur=iQ))]^^[^-^ + Vof[^^-^^ + V^^'icf" = c 

The time the system cost from oscillating to stillness can be bounded by a 
constant. Because of the same reason this method does not break up the lower 
bound of Turing Machine. But for some other things, we tend to pay more 
attention to it. Some relevant issues will be discussed in Sec-V. 

3.6.3 Graph Isomorphism, Graph Spectrum and Oscillators 

In this part of the section, wc shall talk about a complex example in detail. We 
do not mean to show that the method we designed here is superior to all of the 
TMs constructed by the people of the same aim. We just want to demonstrate 
a new style of computation. 

3.6.3.1 Spectrum of Graph 

Suppose X = {V, E) is a graph, A is it's adjacent matrix. We say JaW is the 
characteristic polynomial of X, also denoted by /x(A). (Ai, . . . , A„), the whole 
root of /(A), is called the spectrum of graph X. 

Actually two different adjacent matrices may represent two isomorphic graphs. 
If we alter the permutation of the number of the vertices, A will become P^^AP, 
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where P is the corresponding permutation matrix. However, the characteris- 
tic polynomials of them are the same. Therefore, fx{^) and the spectrum 
spec(X) — (Ai, . . . , A„) are uniquely determined by X. 

For the relationship between spectrum and graph, people conjectured that 
graph can be uniquely determined by spectrum, i.e. suppose 

spec(A) — spcc{B), 

can we conclude that 

Unfortunately, the different graphs of the same spectrum were found soon. 

Nonetheless, calculating the spectrum is also important. Because we can 
know a lot of crucial properties, such as the extensionality, rapid mixing time 
of Markov chains on the graph, by the spectrum of the graph. What's more, 
when two graph have same spectrum, and spectrum is never repeating, we have 
a polynomial time algorithm to check whether they are isomorphic. 

1. Input graphs G'iG'2, compute their spectrum, denoted by A1A2. 

2. Compare the spectrums, if Ai ^ A2, then return NOT ISOMORPHIC ;else, 
continue; 

3. Check whether the product of the two similar matrices is a permutation 
matrix, if it is true return ISOMORPHIC, otherwise return NOT ISO- 
MORPHIC; 

Notation: Here by Ai 7^ A2 we mean after sorting their eigenvalue, the two 
series are not identical to each other. And accordingly Gi, G2 should also be 
altered into Gi, G2. But for convenience, we do not differeciatc Gi and Gi. 
Proof: 

If Ai 7^ A2, then Gi ^ G2. So we only consider the case in which Ai = 

A2 = A. 



i.e. Suppose 
then we have 
thus 



Gi = PAP^, G2 = QAQ"" 
P^GiP = A = Q^G2Q 



:>T\'^ r< fr^TjT\ 



Gi^{QP^) G2{QP^ 



by the preconditon,A is never repeating, so PQ is the unique orthganol 
matrices. 

the rest is to show that if Gi, G2 is isomorphic, then QP^ is the permu- 
tation matrix desired. 
In fact, if Gi = G2,then there exists a permutation matrix S s.t. 

Gi ^= o G20 

Since G2 = QAQ^, the formula above means 

Gi = S'^QAQ'^S = (g'^5)^ A (Q^S*) 
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Because of the uniqueness of P, we can conclude that Q^S = P^ , and by 
orthgonahty of Q, we obtain 

S = QP^. 

n 

3.6.3.2 Harmonic Oscillator of multi-freedom 

Suppose s is the number of freedom of the system, qao{(x — 1, 2, . . . , s) is the 
general coordinates when the system is in balance. Without lose of generality, 
we can always assume that g^o is just zero, i.e. qao = 0(a = 1, 2, . . . , s). 

Because we only talk about little vibration, so we only keep several terms in 
the Taylor series of the Lagrangians L of the system about qao . 

The potential energy: 



y-y'-±(S:)'---±±l 



Note that Vq can be omitted. Introduce the notation ka^ 



which is called the strength coefficient. According to the formula ( ^^ ) — 0, 
the second order of the potential energy could be represented as 



99= /O 



^ s s 

Then assume r^ — ri{q) is not relevant to time, i.e. the obligation is constant, 
so the kinetic energy is: 

ITi ^ n s s i~. i~x 

,,, v^ • • ^ V^ V^ V^ "^J '^''"i ■ • 

^ t=l ^»=la=l/3=l *^^" "'^f^ 

Introduce the symbol TOq/3, 

dri dri 

?Ti, 

then the kinetic energy could be represented as 



ruafi = mi3a = y^ mi—— ■ ——, 
"^ '^ dqa dqp 



^ S S 

T ^ l^'^^mati<io.q(3- 

a=l,3=l 

Keep the formula above to second order and since qaqfi is second order m^^ 
should be expanded to zeroth order. In other words map could be looked as 
constants, we just take the value of them when the system is in balanced point. 
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So the Lagrangian could be written as 

^ s s 



Thus the Lagrangian equation is 



d d /I 

dt dqa \ 2 



EE 



mp^qpq^ 



dqa 






Mqpqi 



= 0. 



/3=17=1 



^ 2 E '™"^'^'>' + 2 E ™/3"9/3 + 2 E ^"-r'^-r + 2 E ^Z'" 

V 7=1 /3=1 / V 7 = 1 /3=1 



?0 = 0. 



therefore 



Let 



y^ maiiqij + y^ fcgggg = (a = 1,2,...,; 

/3=1 /3=1 



qp=Ape^' (/3 = l,2,...,s). 
Take it into the former formula, we get the linear equations for Aj^. 

s 

y^^jmapX^ + kap)Ap = (a = l,2,...,s). 

/3=1 

If the equations have non-trivial solutions, then following conditions should be 
hold: 



TOiiA^ + Zcii mi2A^ + fci2 

m2lA^ + /C21 m22A^ + fc22 



TOisA^ + kis 

TO2sA2 + A:2s 



TOsiA^ + ksl ms2X^ + ks2 ■ ■ ■ ITiss^^ + kg 



==0 



This is the equations of times s of A^ , and we can get s A^ , denoted by 

3.6.3.3 The characteristic oscillators for a Graph 

Making use of the conclusions above, we construct a specific oscillators for any 
given connected graph. 

Denote the vertices of graph by numbers 1 ^ n, according to any order. 
The mass of a vertex is set Ig. Connect the vertex 1 and n to ends by strings 
whose k is zero by that direction. For the rest, we connect them according to 
the adjacent matrix, i.e. if Aij ~ l(Note that Aij ~ Aji), connect vertex i 
and j by a string whose fc == 1. Let's study the motion of the system: First, 
if two vertices is not connected by string, we have ka/s = k^a = 0. Second, 
the vibration is little, so string is not an obligation. And we take the general 
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coordinates as the usual displacement vectors, so map — m^a — ^ap, where Sap 
is the well known Kronecker notation. 

At last, we obtain the determinant as follows, which is the characteristic 
polynomial of our system. 



A2 + di -Ai2 
-A21 A2 + d2 



'Asi 



^A 



s2 



-Ais 

-A2s 

X^ + d, 



It has been proofed that A^ < 0. So let —A = A^ ,we can see that the determi- 
nant above actually compute the spectrum of A' which is converted from A by 
adding inulti-loops(the number of degrees). If a vertex is in the characteristic 
position, it will take part in the vibrations of all frequencies, if no one is in the 
characteristic position, then they just vibrate with respective frequency. In both 
cases, we'll measure the frequency and differentiate them by means of FFT, so 
as to get the spectrum of A' . 

Apparently, adding multi-loops is not harmful to the decision of whether A 
and B are isomorphic, for if A' 7^ B' , then A ^ _B If A ic _B, then A' ic B', 
which will also be checked by the oscillating system. 

3.6.3.4 Comments 

Suppose now the problem we want to solve is: try to find the riQ— th(no is a 
constant) value of such matrices. It is easy to see that the period of the system 
satisfies 

-^<T< ^(or ^/o < / < V^fo), 
V"/o Jo V"- 

where /o is the eigenfrequency of a single string, considering the minimum case 
occurs when the corresponding graph is totally parallel connected(two vertices 
with n— multiple edges between them), while the maximum case occurs when it 
is just a chain. 

The total steps of sampling should be 



N 



2B. = 2(^/,-^)(c„.f) 



where the constant c„„ is related to the required precise uq. So we have: 
0{N) = 0{n) and the corresponding complexity for fast fourier transform 
should be 0(n\nn). 

So we can say that the time complexity of this method should be 



OK) + O ( -^ j + O(nlnn) + 0(n) = O(n^) 

where the left 0{n^) is the cost of constructing the system and 0{n\nn) is 
the complexity of FFT. On the other hand, if we use the well-known algorithm 
called QR-method to get the answer, it will cost such Turing machine 0{n^) 
steps. However, we know little about wether our method is superior to any most 
efficient Turing machines. 
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However, such analysis may not be enough. We cannot obtain a rehable 
result unless some specifications for the materials are considered. The procedure 
of constructing the system according to an arbitrary graph is actually quite 
tricky. For this system need n oscillators of the same mass and different length. 
To achieve this, we have to assume that there exists a kind of 'ideal' material 
which, at least for a large enough range, can be stretched to a desired length 
easily(in 0(n) time and totally 0(n^)) and, at the same time, keep rigid. 

3.6.4 Steiner Tree Problem 

Steiner Tree Problem is a problem in combinatorics. The general version of 
Steiner Tree Problem is NP-complete, which implies that this problem is unlikely 
be solved in polynomial time. 

This problem is similar to the Minimal Spanning Tree Problem in metric 
space. The difference is that Steiner Tree Problem allow people to add new 
points v'{v' ^ V) and new edges e'(e' ^ E) into the original graph G, if neces- 
sary. When \G\ = 3, the new point (in this case, at most one point is needed)is 
called Fermat point. 

At a time, some people became to believe that the experiments of soup 
membrane can be used to solve the Steiner Tree Problem. In fact, when \G\ is 
small, say, less than 5, this method really works. However, when the number of 
vertices is 10 or more, this experiment just cannot give the right answer. One 
can attribute the failure to different reasons and derive various explanations. 
Of all these potential explanations, the one which states that 'it is just the 
errors during the experiment cause the failure' made many people conjecture 
faithfully that classical mechanics can be used to solve NP-complete Problems 
in polynomial time(So they try to proof P=NP). 

In fact, the foundation of the experiment is the well-known property that the 
membrane will stay at a stationary state, where the surface it produces will be 
just the minimal surface. Unfortunately, this theory has nothing to do with the 
fact that the membrane can arrive at the stationary state fast. What's more, no 
one can proof the soundness of such property under the framework of classical 
mechanics. 

3.6.5 DNA Computation 

In 1994, Adleman used a probabilistic DNA algorithm to solve HP problem 
(Hamilton Path Problem). HP problem is NP-complete, which implies it is 
difhcult to find a polynomial algorithm to solve it|14|. 

In order to understand Adleman's method, the following knowledge seems 
necessary. 

(1) DNA contains chains consisted by four types of nucleotides, denoted by 
A, C, G and T. 

(2) These nucleotides forms complementary couples, i.e. A and T are comple- 
mentary, C and G are complementary. If the corresponding positions of 
two DNA chains are complementary, they will patch up as the twin-helix 
structure. 

(3) PCR, which proposed by Kary Mullis, is method to reproduce the specific 
chain we need. 
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(4) There is a machine caUed 'sequencer' which can be used to read out the 
series of a DNA chain. 

Adleman's Algorithm contains five procedures(Suppose \G\ — n): 

(1) Randomly produce the paths in the Graph, encoded by DNA chains. 

(2) Keep only those paths which begin with Vin and end with Vout- 

(3) Keep only the paths whose length is n 

(4) Keep only those paths which enter all vertices in G at least once. 

(5) If any paths remain, return 'True', else return 'False'. 

Note that the first step of Adleman's Algorithm, which is usually thought to 
be work as an initialization operator V, is not polynomial with respect to the 
resource mass 9Jl and space 6 at least. Considering asymptotically we can 
only sequentially get the mass the algorithm need, so actually 0{n\) mass can 
cause 0{nl) time T. As a matter of fact the other steps of this algorithm, 
which require exponentially molecules fully blend by polynomially increasing 
contacting facades, also cost a lot of resource T. 

It is not very hard to appreciate the conclusion that we can obtain great 
power of computation suppose we are provided with corresponding quantity of 
mass, and do not take the cost of preparing such equipment at all. For one 
thing, let's consider the following ideal model. 

Suppose we have enough universal Turing Machines, each of them are de- 
noted by their footnotes. What's more, by some altering in the definition, these 
UTMs have the ability to transmit their results to others. And the condition 
of two UTMs Ui,Uj{i y^ j) could communicate to each other is that they are 
adjacent to each other, denoted by Adj{Ui, Uj). 

So the computational model constructed following, called 'Turing Tree', can 
exponentially speed up the computation of any NP-complete problems. 

Definition 3.14 (Turing Tree) Suppose we have infinite many UTMs, each of 
them denoted by unique footnotes, and 

Adj{U^, Uj)-^j = 2i + l\/i = 2j + l\/j = 2i + 2\/i = 2j + 2, 

then we call this Turing Tree. 

It is easy to see that the following relation holds: 

Adj{Uo,Ui),Adj{Uo,U2) 
AdJiUi,U3), Adj{UuUi),Adj{U2, U5),Adj{U2, Ue) 



For example, a TSP problem can be solved as following: 

a The Observer input the weighted complete graph G to the Uq, Uq decode 
to a permutation and compute the sum of the weight, and then transmit 
G and flag i^ = to Ui,U2. 
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b For index zAfter Ui get F = and G, it check whether i < [logj nl~\ ,if 
the answer if 'yes' then decode i to a permutation and get the sum, and 
transmit G and i^ = to U2i+i and C/2i+2;else check whether i — [log2 nl~\ , 
if it is true, decode i to a permutation and get the sum, then submit the 
weight sum to the C/[i/2]-i- Else, do nothing. 

c For index i, after Ui get _F = 1 and two sum(come from 1/21+1,^21+2), 
if it's index is not zero, then submit minis'^, 5*2;+!, 52^+2} and F = I to 
C/rj/2i-i- Else return minjS'o, Si, S2} and write it on to the tape. 

It is easy to check that the subprocedure of the algorithm which is used to 
decode a natural number to a permutation is polynomial. So the cost of time 
the Turing Tree consumed should be 0(2 log2 n!) < 0{2n log2 n) (Including once 
sharing the task and once championship for the minimum). So it is the time to 
answer how can we 'easily' construct a big enough Turing Tree. 

3.7 Preliminary Discussion of the classic theory of Com- 
putation 

3.7.1 Turing computable is physical computable 

The topic about the existence of a theoretical physical system which can provide 
an implementation of universal Turing Machine has been studied by many schol- 
ars. In addition to the current implementation of computers, scholars have con- 
structed many other wonderful designs on various axiom systems of physics(e.g. 
Classical Mechanics, Quantum Mechanics). 

Of course, the results above only imply that it is the ideal mathematic model 
for a family of physical phenomenons can be look as equivalent to UTM in 
terms of computability. After all, we cannot know for sure that some theory of 
physics is completely correct. Because of this, when we talk about the ability 
of computation for a certain family of physical system, we always assume either 
of the two preconditions: 

• The ideal mathematic model of some branch of physics is believed to be 
absolutely right. 

• At least in a very large scale, the theory works. 

3.7.2 PLATO Machine 

For several decades after the Church- Turing Thesis was proposed, people failed 
to find a counter-example. This kind of counter-example, if they really exist, 
should satisfies the property that most people think they can be effectively 
computed in principle, and no Turing machine can compute them. 

However, many physicists tend to make efforts in another direction, that is, 
they want to find a family of processes in nature, whose functional expression 
may not be intuitively computable, nor Turing Computable, but actually it can 
be used to 'compute' a nonrecursive function by measurement. 

Suppose the problem we attempt to deal with now may cost infinite many 
steps for some computational model(e.g. Turing Machine), does it necessarily 
mean that we have to wait infinitely long time to get the results? This is 
not always the case, PLATO Machine, which was proposed by H.WeyliTj, is 
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just a counter-example. Though it is named after Plato, the designer's main 
inspiration comes from one of Zeno's Paradoxes. 

Specifically, PLATO Machines use (1/2)" seconds to execute the n— th step. 
For instance, suppose the decision problem we want to solve is 3nP{n), where P 
is a predicate and P{x) is used to describe some properties of x. Then PLATO 
machine P will check whether P(l) — 1 holds in 1/2 seconds, and check whether 
P(2) = 1 in 1/4 seconds,. . . , and check whether P{n) holds in 2~" seconds, and 
so on. It is easy to conclude that if P find an answer, it will return the answer 
in one second, otherwise it will return false after a second. Considering the sum 
of geometric series, the proof is trivial. So the upperbound of the time for P to 
solve any question is 

Apparently, if P does exist, its power is extraordinarily great, for it can even 
solve Turing's Halting Problem in one second. 

So far we have seen two idea to implement the PLATO machine P. However, 
unfortunately, neither of them are successful. The first one is to construct the 
machine according to the definitions of H.Weyl. Apparently, it is difficult, for 
people do not believe that time is infinitely divisible. The second one is to make 
use of the theory of general relativity. However, the computing system will also 
exhaust the resource of the universe which make the observer cannot get the 
answer. 

3.7.3 Recursive function whose derivative is not recursive 

April 1970, J.Myhill published his astonishing result[6]: There exists a recur- 
sive function, whose derivative is not recursive. In order to understand the 
principles of the construction, knowing the following fact about the recursive 
functions (whose domain is M) should be helpful. 

Theorem 3.15 Suppose f is a real-valued function, {/„} is a series of recur- 
sive functions, if there exists a recursive function e : N — > N s.i. Va; G //s > 
e(n) \fk{x) — f{x)\ < ^, then f is recursively computable. 

J.Myhill's idea is to build a non-trivial structure(slope or bump) in the 
neighborhood of 2~" in interval [0, 1], where n S £/, and jz/ is a recursively 
enumerable, nonrecursive set. Otherwise f{x) = 0. However, in order to make 
the function computable, the scale of the structure should shrink as the n is 
enumerated recursively, or rather, should be smaller than the bound in the the- 
orem above. As a result, the derivative of the function is intuitively hard to 
compute, and on the other hand we can proof that it is not recursive, because if 
we could compute it we can use the result to decide whether [a;] is an element 
of ^ generally, contradicting the nonrecursiveness of £/. 

Specifically, suppose 

AM-/ a;(x2-l)2, if-l<x<l; 
^^ > - \ 0, if |a;| > 1. 

It is easy to verify that 9{-l) = 9(0) = 0(1) = 09'{-l) = e'{l) = 06''(0) = 
land 9min = 0{—i/V^) = —^Omax = 0{-\-l\^) = -\-X. We call 9 a bump of 
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length 2 and height A. Then the function 0q^(x) = {(3/\)d{x/a) satisfies the 
following conditions: 

0M-a) = 0,^(0) = 0a,p{a) = 0, 9'^pi-a) = ^^(a) - 0, ^^(0) = 0/Xa, 

-I3<0ap{x) < 13 {a<x<a.) 

For each n S jz/we shall construct abump: Oa„p^ at 2~" i.e. 
if n e i2/, (5 e [-a„, +a„], /(2-" + 5) = 0c"pA5), otherwise f{x) = 0. To 
make / well-defined, parameters a„,/3„,ri G ^^ is defined as 

where n ~ h{k) and /i is a function enumerating jz/ without repetitions (It is 
easy to proof that if there exists a recursive function enumerating ^, then there 
exists such function with no repetitions). 

For physicists, does J.Myhill's results imply that if an object move under the 
condition that the displacement and the time satisfies the following relations 

'"'^ ' ■'^ ' \ 0, o.w. 

The speed v = r'{t) will be a physical quantity which is not computable? 

3.7.4 Physical States which is not computable 

Pour-El et al published their results in 1997: for a differential equation, one can 
design a specific initial state to make the solution after t {t could be take some 
computable value)seconds is nowhere computable[8]. 
Consider the IVP of the following wave equation: 

u{x, y, z, 0) = f{x, y, z), ^{x, y, z, 0) = . 

where {x,y,z) eM?, i e [ 0, +oo ) for aU / e ^^^ this IVP has a form of 
solution known as Kirchhoff's formula: 

w(^,t)= // [f{~^ + tn)+t\7f{~i + tlt)-lt]d<7{lt) 

The conclusion Pour-El get is the following theorem: 

Theorem 3.16 For all compact set Z? C K'^ x [ 0, oo ), there exists a com- 
putable function f{x,y,z) € ^^, s.t. the corresponding solution u{x,t) is not 
computable in the neighborhood of any point in D. 

Pour-El et al construct the initial value through the uncomputable real number 

Apparently, one can conclude that in this wave equation, the initial state is 
computable but the state m(0,0,0, 1) is a state which can not be compute. 
For us, can we safely conclude that 

{Turing Computable} C {Physical Computable} 

but 



{Turing Computable} 7^ {Physical Computable} 
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? 



3.7.5 A few Comments 

In the above scenario, the use of (actual) infinity is their common theme. They 
ask the system to run for infinite steps or just encode the solutions into real 
numbers. It is easy to find out that adding either of these two assumes into a 
physical system will make the original system extraordinarily powerful. 

For example, we can throw a particle onto a plane [0, 1] x [0, 1] at ran- 
dom(obey the uniform distribution), then we can proof that with high probabil- 
ity, the x-coordinate(or y-coordinate) of the center of the particle will indicate 
a non-recursive real number. In fact, in cell [0, 1] x [0, 1], the Lebesgue measure- 
ment for the recursive real numbers is 0, while the rest is 1, i.e. 

m([0, 1] X [0, 1] n Rr) = 0, m([0, 1] x [0, 1] n M^) = 1 

This is geometric probability and consider the uniform distribution, the proba- 
bility of the either event of the two are just their measurement. Therefore we 
can look the a;— coordinate as a function with respect to the digits. According 
to Beggs et al, a theoretical machine called SME may help us to get the value 
of the position coordinates [20 ] pi] | 22 ) . 

However, does the strict plane really exist in the physical world? We just do 
not know. 

We propose the some levels for / which is not computable. Suppose o is an 
operator, 6 is a physical implementation of o and operators A , V always exist. 
For f (z F where all elements in F are non-recursive functions, we have the level 
of existence as follows. 

• Existence-I 

3p(Aopo V = /) 

• Existence-I* 

3p3V3F((/ e F) A (Pr{A o p o V = /'|/' G F} > 0)) 

• Existence-II 

3A3p{AopoV = /) 

• Existence-II* 

3p3V(Aopo V = /) 

• Existence-Ill 

3p3A3V(Aopo V = /) 

According to the levels we proposed above, assume the space is continuous, we 
can find out that P G Existence-II, J.Myhill's function / G Existence-I, Pour- 
El's construction g Existence-II, our example x G Existence-I*. Apparently, 
we wish to get the examples in Existence-Ill. 

4 Physical Resource Complexity for Quantum 
Computation 

4.1 Physical Resource Complexity for Quantum Compu- 
tation 

For general quantum computation, we only need to explain the definition of the 
physical state set and the required evolution operators. More over, we only talk 
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about Monte Carlo styled quantum algorithms. 

According to von Neumann's four postulates for quantum mechanics, we 
require that the state of any representation should be vectors in Hilbert space, i.e. 

fi CH 

and the evolution operators should be unitary, i.e. 

H e U{n) 

Without loss of generality, we can assume that the measurement operators 
is projection operators(POVM could be substituted by projection operators 
through adding more auxiliary qubits) 

4.1.1 The RCEF for Quantum Computation 

The resource cost by a computation is 

9a^(Ao-HoV) = (T(Ao-HoV), 
e(Ao-HoV), 
€(Ao-HoV), 
0(Ao-HoV)) 

Our definition here is special a case of the one in Sec-III. Suppose our discussion 
is restricted to QCM, i.e. we have finite kinds of universal quantum operators, 
then the number of gates used and the depth of the whole circuit will be the 
main parameter which should be took into account. It is easy to find out that 
this definition is similar to that of quantum circuit model. One of the difference 
between them is that we will also take the cost of design(usually this costs time) 
of a new circuit into account. Though in most cases, this will not cause great 
difference from the result given by QCM, however, we don't think we can safely 
ignore the potential exceptions just because it is usually easy to expand the 
scale of some circuits. 

So far, people always assume that qubit is relatively easy to prepared. At 
least in the asymptotic sense, no matter how difficult to prepared a quantum 
bit, the cost should be bounded by a constant. We will also do this. 

4.1.2 Deutsch-Josza Algorithm 

Deutsch-Josza algorithm is one of the most successful algorithms in the early 
years. The corresponding problem of the algorithm is: consider two sets of 
functions: 

A: L\^ : {0, . . . , 2" - 1} ^ {0, 1}, Vx(^(x) = 0) 

B: L|^:{0,...,2"-1}^{0,1}, |{x|^(x) - 0}| - |{x|^(x) = 1} 

Apparently we have AO B = (d, now suppose f E AU B and there is an oracle 
to compute /. We arc required to decide whether f E A or not. It is no doubt 
that people wish to reduce the times of query the oracle as much as possible. 
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Note that the cost of implement the oracle is not taken into account, because 
we assume we have implemented it. 

The algorithm needs a trivial input ipo = jO)**"!!), and used the gate iJ'*"(g) 
H onto the stateVJo a nd get V'l , i-C- 

V-i = (^^"(^if) (|0)®"|i)) 

Note that H = -^({\0) + |1»(0| + (|0) - |1))(1|). By induction we have 

a;, J 

where i.e. x • y = Xi A j/i. So we get: 



^,y 



H^^^H] (|0)^"|1)) 

( S-(-^\^■v\x){y\]\Q)®-■ 






|o)-|i) 

V2 



(by orthogonality) 



^E(-l)°l^> 



|o)-|i) 
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2"-l 

1^ E N> 



|o)-|i) 



Now use the oracle Uf : |x, y) — > |x, y © /(a^)) onto the state -i/"! to get V'2 



^^2 = t// ^^ E 



|0)-|1) 
V2 



= T^ e'|:^>(/(x)®M) 



2;=0 
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1 '"^\_1)/(.)|^) 



At last we use H^' 

1p3 = 



x=0 

I onto V'2 to getV's: 



|o)-|i) 

72 



= rff«5 



'^)(:7k'E'(-iK(^^l^) 



x=0 



|0)-|1) 
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(-1)" 



= EE 

z X 



^ + /(x) 



k> 



|0)-|1) 
72 



The observer is supposed to check the first n qubits, note that the amplitude of 
|0)®" is J2xi-^y^''^/^"- If / G ^: fix) is constant and the amplitude of |0)^" 
is +1 or —1. So the amplitude of another cases should be zero and the observer 
wiU get jO)*^". On the other hand, ii f e B, the amplitude of jO)*^" will be zero. 
So the observer will always get a non-zero vector. 

In our opinion, the procedure could be written as follows. 

V = Initialize the state|0)'^"(g) 1 1) 

Generate the whole circuit 
Hi = if ^" (g) H 
H2 ^ Uf 
Hy, = H'^"' ® / 

A ^ El^^)(^.l 
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Let -H = -Ha o ^2 o Hi 

(A o -H o V) (|a;)) = (a o % o Hz o Hi o v) (|a;)) = P(/ G B) 

Though Deutsch-Jozsa Algorithm is great, someone still think it is not very 
useful. In addition to the fact that the problem they studied is not very impor- 
tant, there does exist an efficient classical probabilistic algorithm to solve the 
problem with high probability. 

4.1.3 Shor's Algorithm 

Shor's Algorithms for prime factorization and discrete logarithms [TOl [H] is so 
far the most exciting quantum algorithms. The appearance of Shor's Algorithms 
is the greatest challenge to strong Church Turing Thesis. 

Shor's Algorithms depends on a technique of so called "quantum Fourier 
Transform" . But of course QFT is not enough. Shor's Algorithm is totally 
non-trivial and marvelous, and few people can produce any algorithms like that 
easily. 

In order to understand Shor's Algorithm, it may be enough to gain a clear 
idea of quantum ordering algorithm. This is the only subprogram in the Shor's 
Algorithm which has to be implemented by quantum computers so far, and it 
is really the most important subprogram. 

First, note that 

r-l 

y ^exp{—2nisk/r) ~ rSko 

s=0 

and deffiie lu,) as follows 



r-l 

^ ^ y e-""^"- /"^ a;*^ modiV 

^ fc'=0 



According to the fact above, we can get 

r-l 



/r 

s=0 



In fact 

r — l _, r— 1 



1 ^^2.,.fc/,.|^^^^ 



„*: 



mod N 



1 ^g2...feA|^^^ ^J-J2 ( e^--^/^^ Y. e-'"^"'/-^ x'^'modN) 



s=o ^ s=o \ ^ A;'=0 

r-l 

k' 



modiV) 



\ V^ I 2-nisk/r V^ -27risfc'/r 

s=0 \ fe'=0 ■ / 

;Ei:=xp(?=J(i^)|."„.od») 



fc'=Os=0 
1 '■"^ 

- ^ r5kk' x'' modiV ) 



r 

k'=0 



X mod N 
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In particular, when fc = 0, we have 

r-l 



\(g,L 



where L = [log(A^)] . 

Suppose Ux,N satisfies Ux.N\y) = |a;y(mod7V)). Considering Z|J^ and the 
fact that the permutation on orthnornial basis can be represented as a unitary 
operator, one can know for sure that Ux.n is unitary. What's more Ug is a 
eigenvector of Ux,n , the corresponding eigenvalue is e~^ since 



Ux,n\Us) = 



1 ''"^ 



fc=0 



c'^+^modN) = 



Reverse the results above, we get the first half of the quantum ordering 
Algorithm, which complete the following task: 



\«iL 



T^El 



u'J 



Modular exponentiation 




x'^mod N 



where the state ^ 



~^ X]s=o e ''■jUs) is the one we desire. Apparently the eigen- 
value contains the information of r. So as to extract the information, we need a 
sub-progress named " quantum phase estimation" which based on inverse quan- 
tum fourier transformation. One can verify that if t is large enough, such like 
t = 2L + 1 + [log (2 -I- ^)] , for each s S {0, . . . , r - 1}, we wiU obtain the 
estimation oi ip Ri s/r accurate to 2L + 1 bits with probability at least (1 — s)/r. 
Through the continued fractions algorithm, we will get r with high probabil- 
ity (According to PNT). 

In our opinion, the procedure above could be written as: 



V = 

Hi = 

-Hi = 

A = 



Initialize the statclO)*^* <Si 11)*^ 
Generate the whole circuit 

CUx,N 

FT^ (g) /«'-^ 
CF ® I®^ 

Y.\P^){P^\ 



It is easy to check that except the "Ha, all operators cost polynomial time with 
respect to log N . The complexity of operator modular exponentiation and con- 
tinued fraction are both 0{L^), which are two most time-consuming subproce- 
dure of the whole algorithm except the "Ha (inverse quantum fourier transform). 
Note that "Ha is indeed not an operator which could be implemented by 
polynomial universal gates. Consider a family of gates used in "Ha which is 
usually noted by Rk{k € {2, . . . , L}) 



Rk 



1 

e 





27ri/2'' 
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In other words, the original Shor's Algorithm is not a algorithm with super- 
polynomial acceleration. In order to overcome this, Coppersmith created a new 
algorithm called the AFFT(Approximate Fast Fourier transform) [T^j which can 
substitute for the procedure QFT. 

4.1.4 Grover Algorithm 

Quantum Search Algorithm [T5], also known as Grover's Algorithm, is another 
quite successful quantum algorithm. Though this algorithm is not faster than 
the fastest classical search algorithms super-polynomially, one can proof it is 
the fastest one considering quantum mechanics. Therefore, the complexity of 
the algorithm is the complexity of the problem it deals with. 

The crucial subroutine of Grover's Algorithm is the Grover iteration, often 
denoted by G: 

• Apply Oracle O : \x)\~) -^ (-l)^(^)|a:)|-) 

• Apply Hadamard GatesiiJ**" 

• Perform a conditional phase shift(2|0) (0| — /) on the computer, with every 
non-zero bases receiving a phase shift of —1. 

• Perform Hadamard transformation H®'^ . 

Note that iJ^"(2|0)(0| - /)i/®" = 2|?/')(V'I - ^ One can proof that Grover 
iteration can be looked as a rotation in the plane spanned vectors which denoted 
the right answers and the wrong answers. 

Let S^ be the sum of all the vectors which indicate a solution to the search 
problem, E^' the rest. Define normalized states: 

1/3) ^ -7M^» 
thus the initial state ji/') — -j^rj2^x=o\^) could be represented as 




The action of Operator O is 0{a\a) + b\fi)) — a|a) — ^1/3), which could be 
looked as perform a reflection in aji— plane. Similarly Operator 21-0) (^| — / 
also performs a reflection in a/3— plane. Thus both two reflections which could 
be looked as a rotation occur in the a/3— plane. Let cos 0/2 = \/{N — M)/N, 
s.t.|-0) = cos0/2|a) + sin6'/2|/3), apply the iteration once makes \ip) become 

G\^) =cosy|a) + sin y 1/3) 
k times use of Grover's Iteration will lead to the following result: 

.'2k + l^. , . nk + i 
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Since \ip) = ^J{N - M)/N\a) + y^M/N\l3), we just need to rotate IV') 
arccos y/M/N radians to the one which is parallel to vector |/3). So repeat- 
ing G ioT R = [ J- ] times will get 1-0) to within an angle 6/2 < 7r/4 of 

1/3). This is a 'good' state, for people only have to repeat the experiment for 
expected constant times to get the solution to the problem(Consider geometric 
probability distribution: E[X] = 1/(1/2) == 2). 

Apparently R < [7r/26'] , suppose M < N/2 then we have f > sin | = J^. 
Thus, we obtain: 

R< 




in other words we need repeat G for i? = 0{y/N/M) times. 

4.2 Quantum Simulation and Quantum Algorithm 

Quantum Lattice Celluar Automata(QLCA) and Quantum Gas Automata(QG- 
A) are two familiar ideal models in the research of quantum simulation^lS . 
Meyer, Boghosian[TSl [THl [H] have obtained their results respectively by using 
these models, that is, they construct some quantum algorithms which demon- 
strate exponentially speedup in such models. For Bohosian, the object they 
tried to simulated is a QGA which obey lattice Boltzman distribution, where 
arbitrary fields can be concerned with. They have proofed that the complex- 
ity of simulation is only related to the dimension of the lattice, but almost 
has nothing to do with the number of the particles. However, the number of 
particle always cause exponentially hardness on a classical computer. In fact, 
Boghosian's results imply that it is almost impossible for a classical computer to 
simulate one evolution step of a quantum system including dozens of particles. 

We've mentioned that it is the difficulty of quantum simulation that makes 
people believe quantum mechanics can provide enormous power of computation 
in the early years. 

Note that in this article we do not care about the hardness of simulations. 
Generally speaking, the hardness of simulation has nothing to do with the one 
of computation. For instance, people may find it difficult to simulate some 
classical celluar automaton according to the given regulations, however once the 
tedious work has been completed there often exists some more simple methods 
to produce the series. A typical example is that the regulations of an automata 
actually cause a circle with a finite period in the series. The same thing can 
happens to quantum simulations too. 

However, it is important to know that there must exists some cases in which 
simulations and computations are equivalent. These extreme cases often ap- 
pears when the length of regulations is near the Kolmogorov complexity (lower 
bound of description) of a series. Still, strictly speaking, at present no one 
can proof that polynomially universal unitary operators really cause exponen- 
tially difficulty in classical computation. To understand this, just consider an 
easy but helpful fact that almost all the problems we want to efficiently solved 
on a quantum computer are in the class BQP, and we have BQPCPSPACE. 
Unfortunately PSPACE=P is not totally impossible. Of course most people 
don't believe this is true, since this would imply that Shor's Algorithms can be 
polynomially simulated on a classical computers. 
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Now we discuss how to extract a corresponding quantum algorithm from 
a method of quantum simulation, which is believed to be exponentially faster 
than any classical one of the same target. 

On a high level, we should do following things: 

• Find a family of experiments of quantum mechanics which can be effi- 
ciently simulated by quantum computers but are believed to be hardly to 
simulate and compute by classical computers 

• Design a 'good' problem about some non-trivial properties of the last state 
of the system, which makes quantum computers able to present the answer 
to the observer quickly. 

Designing the problem is a crucial step. In most cases, though we may have 
quickly obtained the probabilistic distribution very close to the real experiments, 
we can not know the whole information in short time. So first we have to ask 
a question which can be easily verified by any quantum computers containing 
the whole quantum information of the system. 
For example, we can ask a question such like: 

• What the number of the no— th digit of the probability of a certain system 
arriving in 51' (17' C il)? 

The problem of this method is that in high dimensional spaces, it is very 
likely that the probability of the set ft' is exponentially close to zero, which 
actually enables a classical computers to guess zero without running and get 
the right answer in most cases. 

Now we propose our version: Suppose (j) is the wave function of the system 
we've simulated and \(p{X)\'^,X C il is the probability of x appear in X. Try 
to find two subsets A, B C fl s.t. 



3 ^ \HAW 
7 - \HBW 



< 1 



and determine the value of the no— th digit of 0(^4). 

For the systems which (probabilistic) Turing Machine cannot simulate in 
polynomial time, the question above is intuitively hard to answer, though up 
till now no one can proof or disproof it. 

On the other hand, if these systems can be efficiently simulated by quan- 
tum computers, repeating following procedure will ensure us to find the answer 
relatively much faster than any probabilistic Turing Machine of the same aim. 



Definition 4.1 (Vector of normal vectors x) 







Vl/ 



1 





voy 



Definition 4.2 (Procedure Pq) Pqiln pseudo-code): 

while(find the answer) 

{ 
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Mid-cut the space H. by super-plane whose normal vector is Xi 

Suppose the two spaces is fii and J72 

if (the condition is satisfied (verifies by testing)) 

{ 
halt 

} 
else 

{ 

0. = min|0|{r2i,ri2} 

i + + 
} 



5 Conclusions and Future Works 

We formally proposed the theory of physical computation, define the concepts 
of resource and complexity. Several examples, including classic mechanics and 
quantum mechanics, were discussed and analyzed under the framework of phys- 
ical computation. A technique, which is used to converse a method of quantum 
simulation into a quantum algorithm, is discussed. 

This is a exciting field, we believe there is more exciting topic to discussed. 
A very interesting question is: can we find a physical mechanism as the fastest 
implementation of an arbitrary functions? 

In Sec-III, we talk about the question of calculating the centroid of an ob- 
ject. We thought it is the limitation of dimensions (only three dimensions) hide 
the advance of the method we mentioned. We conjecture that this method 
has a excellent counterpart in high dimensional cases. We'll have a try in the 
(quantum)statistics mechanism. 

In Sec-IV, we talked about quantum simulations and how to construct a 
clever problem to induced a quantum algorithm. Actually, we conjecture that 
the problem we construct is a hard one in class #P, for these questions have a 
counting style. However, we are not sure about whether the designed questions 
could be in #P — hard under some specific statistical models. We shall try to 
work on this in the future. 

We've mentioned that we assume that polynomial qubits is polynomially 
hard to prepare. However, it is harder to control the qubits as the number of 
them increase [TS] so far. So one can still conjecture that preparing qubits itself 
is a " complicated computing" , and the results up till now can be explained as 
someone displace the resource consuming procedure, just like DNA Algorithms. 
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